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RETRIEVAL OF DISPERSIVE AND CONVECTIVE TRANSPORT PHENOMENA 
IN FLUIDS USING STATIONARY AND NONSTATIONARY 
TIME DOMAIN ANALYSIS 


SECTION I. BACKGROUND 


A. Introduction 

New algorithms have been developed which permit the separation and consistent 
interpretation of the simultaneously occurring kinematics — convection and dispersion — 
associated with stationary and nonstationary fluid transport. 

The separation of the fluid kinematics of the transport process poses both an 
experimental and an analytical problem. The experimental problem is to retrieve sufficiently 
detailed local information for a field description of the scalar quantities without introducing 
probe interference on the field or the modes of transport. The solution utilized was to 
remotely sense by optical means — that is, using radiometers — a local flow field by 
crossed-beam triangulation. The temporal information thus retrieved contains the signature 
of the space^time variations, of the transmission coefficients (absorption, scattering, and 
emission)[ 1 ] . 

The analytical problem is to extract a description of the transport kinematics in 
terms of the field parameters that are related to the accepted description of turbulence 
flows. A description of the convective component of the transport can be retrieved from the 
temporal histories of the fluctuations by using two radiometers and cross-correlating these 
histories. This cross-correlation results in the average temporal relationship between events 
occurring at two different points along the streamline. This transit time is used in 
conjunction with the transit distance — the distance along the streamline between the points 
of observation — to calculate the convection speed. This experimental and analytical 
procedure is called the crossed-beam correlation technique. 

Fluid transport investigations with the crossed-beam correlation technique have been 
conducted in cold jet flows [2], in tropospheric turbulence flows [3], in glow 
discharges [4] , and in combustion processes [5] . The analysis used to interpret these 
experimental results accounted for only the convective component of the transport.- 
Consequently, an incomplete description was obtained in that it failed to account for the 
dispersive transport component. This analysis also failed to afford a systematic classification 
procedure and failed to establish the probabUity for the transport retrieval occurring along a 
common streamline. 

The description developed here compresses the temporal fluctuations into a wave 
packet; whereas, the presently accepted description uses an infinitely long wave description. 
A wave packet is recognized by its finite coherence length, whose value can be determined 



from correlation analysis. The key to the separation of transport phenomena is the 
recognition that the characteristics of the wave packet obtained from the finite data 
record contain the signature of the component of the transport phenomena as viewed in 
a moving reference frame. The accumulative time averages obtained from the summation 
of piecewise estimated statistical averages will provide results that differ according to the 
sensitivity of the accumulation procedure for the various transport phenomena. Analysis 
of data sequences in terms of the characteristics of the wave packet and the temporal 
averages of the accumulative procedures results in enhanced signatures for the desired 
separation of the simultaneous but different transport phenomena. 

Section II reviews the ordinary correlation techniques used for stationary data and 
includes additional criteria that have been developed for the recognition and processing of 
nonstationary data in a fixed reference frame. Section III illustrates the transformation of 
the covariance curve into a wave packet description in a moving reference frame. When 
the convection signature is suppressed, the covariance curve is then fitted with an 
algorithm that affords a wave packet description of the experimental time histories from 
more than one probe. The parameters obtained for this wave packet afford a systematic 
classification procedure of the process in terms of a bandwidth-frequency ratio. The rate 
of change of the average frequency in the moving frame gives the signature whereby the 
dispersion is calculated. This unique description affords separation of the components 
associated with the transport phenomena. 

In Section IV, this change analysis is applied to data obtained from both 
laboratory and field test experiments, representing stationary and nonstationary 
ensembles, respectively. Since dispersion tends to make any data record nonstationary for 
large scale waves, a discussion of nonstationary data is mandatory. Considering all data 
that have been evaluated, the effectiveness of the proposed methods is then assessed from 
an evaluation of all the results. 


B. Fluid Model 

An analysis of the turbulent transport process which describes the coupling 
between the mean transport and the turbulent fluctuations enables us to demonstrate the 
need for a semiempirical, statistical technique in the investigation of fluid transport. Our 
fluid model (Fig. 1) assumes a flowing turbulent fluid that is both inhomogeneous and 
stationary. The inhomogeneities are primarily a result of velocity and density fluctuations 
in the fluid generated by the shear stresses associated with the transport process. The 
density, p , and the center of mass velocity, v , can be described in terms of a mean 
component, (“) , and a fluctuating component, ( )' , as 


p(t) = p + p(t)’ 


( 1 ) 
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PATH DESCRIBED BY 



USED AS SPACE FIXED 
CONTROL SURFACE 

Figure 1. Space-fixed control surface in a flowing fluid. 


and 


V (t) = V V (t)’ . 


( 2 ) 


The condition of stationarity implies that the mean values are independent of 
time. This description is applicable to multicomponent fluid transport such as that 
associated with moving striations in a glow discharge — the shear layer formed by the 
exhaust of an air jet, and the atmospheric transport phenomena — all of which will be 
considered later as the empirical results of convection and dispersion are evaluated. 

To analytically establish the significance of a description for the fluid transport 
involving both a convective and dispersive component of transport, we must first define a 
space-fixed control surface. The control surface is assumed to lie along the axis of a 
streamline, where the streamline is defined to be the path of the center of mass which .is 
generated by the motion of a small fluid volume. While this discussion will be limited -to 
the motion of the center of mass, this streamline can be defined more generally so that it 
includes any conservative property [6]. The control surface will be the stream tube in the 
flow selected such that the mass average flux across the mantle is zero. This directly 
implies that: 


1 . The mass flux entering one end of the control surface is exiting at the other 

end. 
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2. In an expanding flow, the area of the exit of the tube is greater than the 
area of the entrance. 

In addition, we can assume that within the stream tube, a local isotropy exists. 

The continuity equation for the fixed system is 

■It = div pv . (3) 


This equation, when written for the stream tube in terms of the mean and fluctuating 
components, becomes 


III 


d[p + p(t)'] 

a t 


dV = Iff 

V 



— ^rp + p(t)'] [v. + v(t)'.] dV. 
9 X. ^ 


(4) 


When the time average is taken, the last equation reduces to 


Iff Iff div(p7 + p'7' ) dV, 


(5) 


since the mean component is time-independent and the time average of the fluctuating 
component is zero. This follows since 


P' = ^ I P' (t) dt = 0 
0 

p' V = ^ V / p' (t) dt = 0 , 

0 


( 6 ) 


(7) 


4 



and 


— T _ T_ 

P V = — p / V (t)' dt = 0, (8) 

0 


if we assumed a self-stationary flow. By utilizing Gauss’ theorem, the right side of the 
continuity equation for the mean motion, equation (5), can be rewritten as the stream 
tube surface integral; 


1 


/// 



2 3 

//p7 ds+ //p' 7’ 
s s 



(9) 


Consider now the significance of the terms of this equation. 
Term 1 


Ilf 


9t 


dV = 0 


since the flow is constant. 
Term 2 


( 10 ) 


Jf p V dS . (11) 

®end 

The term is the flux of the time-averaged to volume concentrations which are 
convected with the time-averaged mass velocity ^ in the stream tube. By definition- of 
the stream tube, the mass transport from this term across the mantle is zero. The integral 
over one end of the stream tube represents the flow measured by a probe in a fixed 
reference frame and will be referred to as convection. 

Term 3 


1 ^ 
If T / 

S 0 


dS . 


( 12 ) 


5 



The average longitudinal component of this term is zero, since the backward expansion of 
a small fluid element will be cancelled by the forward expansion of the next small fluid 
element as a result of the time averaging. This is why a probe, which is a space-fixed 
observer, in the flow cannot measure this fluctuating longitudinal transport. The 
transverse component is the dispersive transport which results in the lateral expansion of 
the stream tube. This term represents the covariance between the convecting and 
convected properties as measured simultaneously at the same spot by both a space-fixed 
observer and an observer moving with the convected fluid. 

This dispersive term in the continuity equation for the mean motion is the term 
that has remained analytically unsolved [7]. The new approach introduced in this report 
to solve the problem is the employment of empirical information retrieved by a remote 
optical detection system which monitors the density fluctuations, and it is coupled with 
new statistical correlation techniques to extract the dispersive information. When the 
cross-covariance curves are obtained for two path lengths along a common streamline, the 
change in the wave packet represents the signature for the mean longitudinal fluctuations 
of the mass transport as seen by the moving observer. Since we are assuming local 
isotropy, this mass transport is the dispersion. To verify that the measurements are taken 
from probes on a common streamline, it is necessary to calculate the transport 
probability. 

Before discussing this probability calculation, it is necessary to define our basic 
statistical techniques. To begin, we wilt examine the single-probe averaging techniques, 
which are fundamental to this information retrieval analysis. 


C. Fundamental Averaging Techniques 

The averaging techniques are used to establish probe moments. The single-probe 
averages used are the mean and variance, and the two-probe average used is the 
covariance curve. Although N-probe averages do exist, they are beyond the scope of this 
investigation. Our immediate attention will be focused on an evaluation of the first and 
second moments for single-probe averages. 

Assume that the behavior of a physical parameter in the fluid model can be 
expressed as the time history of the amplitude fluctuation, x(t) , of this parameter. This 
terriporal history is thus a function of the running parameter t (Fig. 2), the relative time. 
The range of the relative time is from zero to T. The interval between discrete samples is 
giv,en by the sampling interval. At , such that 


T = N At, 


(13) 


where the integer N is the total number of discrete samples in the data record. For 
easier analysis and manipulation, the total data record can be subdivided into equal time 
segments, AT , called pieces, such that 
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CONTINUOUS NOTATION- x(f) 




•al 


■at. 


-at. 


at. 


EXAMPLE OF T, 


m 


ACCUMULATIVE 2 INTEGRATION TIME 


,3 -p .. ,4- 

WHERE: AT, = AT, = AT, = AT, = AT, 

1 Z 9 4 I 


■T„ = T 


TOTAL INTEGRATION TIME 
Figure 2. The notation of the time series. 


T = T = M AT, 
M 


(14) 


where the integer M is the total number of pieces available in the data record. This 
piecewise subdivision of the total data record not only affords more computational 
flexibility, but this piecewise subdivision also affords an opportunity to apply statistical 
techniques for information enhancement and analysis. In fact, this piecewise analysis is 
the cornerstone in the evaluation of nonstationary data records. Each piece of data 
record contains n discrete samples, such that 


N = M n 


(15) 


or 


AT - nAt. 


(16) 


The accumulative results of a summation of a number of pieces of the data are signified 
by the subscripted lower case m . For example, Tj^ means the interval time from ,the 

start of the data record to the end of the mth piece, or Nj^ is the number of data 

points from zero to the end of the mth^piece of data. The accumulative result will be 
precisely contrasted to the piecewise result, which is given by the subscript i . 

The first two moments of the single probe averages are schematically shown in 
Figure 3. The first moment xj is 


7 



MEAN AND VARIANCE 



X. 

1 


1 

AT 


i AT 

f X. (t) dt 

(i-l)AT ^ 


(17) 


which is known as the mean value of the ith piece of data record. The integral is used 
here as a shorthand notation for the summation when the summation represents a 
continuous sum. In electrical terms, xj is the dc component of the time series, 

and [Xj(t) - xj] is the ac component of the time series. 

The second moment xj^ is 

_ i AT 

= f ^i(t)'dt , - (18) 

' (i-l)AT 


which is the mean-squared value of the ith piece of data record. 


The mean-squared deviation about the mean, (Ax)j^ , is 



1 

AT 


iAT 

/ 

(i-l)AT 



x.]2 dt 


(19) 


commonly known as the variance. The square root of the variance is the standard 
deviation (or the rms value of the ac-coupled signal). This gives a measure of the 
dispersion of the data about the mean. 
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We now introduce some subtleties in notation to show how the mathematical 
manipulation is carried out between the different pieces of the data record. If we use a 
direct summation of data points, then the double bar notation is used; e.g., 


■~T 1 

= — 

m T 


-L 

/ x(t)^dt. 


( 20 ) 


If the results of the individual pieces are summed, the notation is 


<X“=> 


m 


m 


m 



( 21 ) 


Usual methods of time series analysis assume that the accumulative average should 
approach a limiting value, a condition referred to as self-stationarity. Section II introduces 
a two-probe test of stationarity which is sensitive to the signatures associated with 
convection and dispersion rather than with the level of activity. 
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SECTION II. COVARIANCE OF SEPARATED PROBES 


A. Introduction 

The one-probe averages just discussed will now be extended to a two-probe 
average known as the covariance curve. It will be demonstrated how the coherence length 
and the transit times are retrieved from the covariance delineations. In addition, 
techniques to establish the confidence level of the covariance curve along with the period 
of stationarity are discussed. These concepts are then extended so that we can obtain the 
probability of retrieving the transport phenomena for a common streamline. 


B. Accumulative Covariance Technique 

‘ Correlation techniques are employed to compare the measurements from two 
probes which are separated in space and in time. In this technique, the average product at 
a time delay between two time series results from one time series being advanced or 
delayed in starting time relative to the other time series. Hence, the correlation curve 
gives a comparative measure of the degree of similarity as a function of time delay [8,9]. 

To accomplish the above product integration, the total data record is subdivided 
into pieces in the manner discussed in Section I (Fig. 2). Each piece of the data record 
affords a piecewise product mean value as a function of time delay. These values are 
usually normalized with the geometric mean of the mean-square value. The result of such 
normalization is called the piecewise correlation coefficient. We have 



T 1 


AT 


x/x? 


2 2 
i 


AT 

/ x (t-r-T) y (t) dt 


( 22 ) 


where the presuperscripts x and y of the correlation coefficient signify that this is the 
cross correlation of the ith piece of the x and y data records. The presubscript t refers 
to the functional dependence on time delay, t , and the postsubscript i refers to the ith 
piece of the data record whose leng^ is AT . The R signifies the correlation coefficient 
over this ith interval. The terms xj^ and yj^ are the mean-square values of the ith 
-interval, equation (18). 

i 

When a data record is correlated with itself, the correlation is called an 

autocorrelation. The piecewise autocorrelation coefficient, , for the ith piece of 
the x data record at the time delay t is 
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1 


(23) 


^R. 

T 1 


AT 




At 

I 

0 


x^(t 


t ) 


x.(t) 


dt, 


which has a maximum value of one at zero time delay. That is, when r = 0 , there is an 
exact match between the data record and itself. The unnormalized autocorrelation 
coefficient is the mean-squared value at t - d . 

The correlation function [equation (22)] has an inherent limitation. When there is 
a large mean value associated with either (or both) the x and y data records, then the 
small fluctuations are swamped out by this large constant mean. When the swamping 
effect exceeds the dynamic range of the systems, the meaningful intelligence cannot be 
retrieved. 

This swamping effect associated with the correlation coefficient can be 
surmounted by computing the correlation coefficient about the mean of the data records. 
This form of correlation coefficient is referred to as a piecewise cross-covariance 

coefficient, , and is defined as 


xy— 

r. 

T 1 


AT 


AT 7 (Ax) f (Ay) 2 0 


/ [x.(t-HT) - X.] [y.(t) - y.] 


dt 


1 ' “ ' 1 


(24) 


where the autocovariance coefficient is 


AT 


XX— 

r. = 

T 1 


/ [X.(t+T) - x][x(t) - X] 

f. \9 y 1 11 1 


dt 


(25) 


AT(Ax)j 0 


The unnormalized autocovariance coefficient at r = 0 is the variance. 

The traditional time average is a straight time integration over the entire data 
record [8]. For self-stationary data, this straight time integration is more meaningful 
since the average is approaching its limiting value as more information is processed. 
However, very often the boundary conditions for the experiment cannot be controlled by 
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the investigator, and the time integration becomes ambiguous because of the changes in 
boundary conditions which can prohibit the averages from reaching asymptotic values 
with an increase in integration time. Accumulative covariance techniques, where the data 
set is subdivided into segments, provide an opportunity to recognize and account for 
these changes in boundary conditions. 

Sometimes a short segment of the data completely dominates the data record and 
prevents an intelligent evaluation of the experimental results. Such domination can be 
reduced by a simple quantization technique where each segment of the data record is 
normalized relative to itself [ 10,1 1 ] . This quantization affords only an estimate for the 
accumulative covariance. The choice of the mean value and the normalization factor 
provides the different accumulative covariance estimates. This discussion uses the 
following definitions for the accumulative covariance techniques: 

1. Sequential covariance where the means and the variances refer to the time 
integration over all the preceding record; i.e.. 


T m 


mAT 


m 

Z - 
J 


AT 

/ 


[x.(t - t) 


X ] [y.(t) - y ] dt 

m ' m 


(Axf (Ayf 


(26) 


where m is the accumulative total of pieces that have been processed. 

2. Quantized covariance where the means and variances are recomputed for 
each segment of the data record before accumulation; i.e.. 



m 


m 


Z 


i=l 


xy- 

T 1 


(27) 


The piecelength of the data segment of the quantized covariance acts as a low 
frequency filter which suppresses the information associated with frequencies of periods 
greater than the piecelength. This is not true for the sequential covariance, since the time 
averages are not restricted to intervals less than the total accumulative integration time. If 
a specific reference is not made to the accumulative technique used, it implies that both 
techniques afford equivalent results. 
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In Section I (Fig. 2), the interval between samples was defined to the At, where 


At = At. 


(28) 


The maximum time delay, will be one-third the piecelength. ' 

To select the piecelength of a data segment, let us examine the effect of 
piecelength duration on a cosine wave. The resultant autocorrelation curve of a cosine 
wave obtained over the finite interval AT is 

= cos (27 t f t) - E, (29) 


where the truncation error, E , is obtained directly from a table of integrals as 


2 sin^ 27 t£AT sin 27rfT 
sin47rfAT+ 47rfAT 


(30) 


This term could introduce as much as 12.7 percent distortion in the amplitude of 
the correlation curve. We would like to remove this distortion caused by the truncation 
error, because this truncation distortion introduces a distortion in the wave packet 
description. Let us examine the conditions necessary for no, or at least minimum, 
-truncation distortion. When the integration time, AT , is infinite, no distortion exists; but 
the data of interest in this discussion are always of finite duration, rendering this 
condition unattainable. When the integration time is 

AT = ~ , n= 1, 2, 3, ... , (31) 


again no distortion exists. The infinite integration time and the multiple period 
integration times are the limits of integration normally found in statistics texts [8,9]. 

Since random processes normally involve data with a finite bandwidth, the 
multiple period integration time is not directly applicable. The problem is that the 
integration time can be selected to eliminate only the distortion from a single frequency, 
thereby leaving a degree of distortion in the remaining frequencies in the bandwidth. The 
solution is to select a multiple period of the mean energy-bearing frequency to reduce the 
principal distortion. As just pointed out, the truncation error can be reduced by 
increasing the integration time; therefore, the distortion caused by the bandwidth can be 
reduced by selecting a relatively large multiple of the period of the mean energy-bearing 
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frequency. Analysis shows that, when n is eight, the maximum truncation error is less 
than 2 percent. Thus, the duration of the piecelength should be at least four times the 
period associated with the mean energy-bearing frequency to minimize the truncation 
error introduced in narrow and broadband data. 


C. Transit Time and Coherence Time 

The covariance curve contains the signatures for the transit time and the 
coherence time. The transit time is that time required by the transport phenomenon to 
pass between the two points of observation [12]. The coherence time is that period of 
observation for which the transport can be treated as belonging to the same wave 
packet [ 13] . 

To demonstrate the concept of transit time, the elementary model is illustrated in 
Figure 4 [14]. The data collection system is an optical detection system [15] comprising 
two independent units (Fig. 4a). Each unit has a photodetector which receives collimated 
radiation modulated by events occurring at the point where the beam and the streamline 
intersect. These modulations form the temporal history of the intensity fluctuations. To 
illustrate the operation of this system, a test object is moved along the streamline to 
produce the sequential temporal histories of the intensity fluctuations (Fig. 4b). A 
correlation of these two time histories provides a covariance curve that indicates the 
temporal relationship between the x and y data records as the x data recorder is 
advanced or delayed in time (Fig. 4c). By inspection of the covariance curve, we observe 
that the maximum correlation between the temporal histories is 2 seconds — the time 
required for the phenomenon to traverse the distance along the streamline between the 
beams. Utilization of this transit distance accompanied by a knowledge of the transit 
distance between the two beams thus affords us the convection speed of the test object. 

While this example of transit time retrieval is rudimentary, it is the conceptual 
basis of correlation statistics and illustrates the principle applied in complex data records 
containing both dispersive and convective information. A more typical data record that 
can be expected in this work is the initial example given in Figure 2. 

The coherence time is obtained from the covariance delineation. The coherence 
time, which is a direct result of the application of the uncertainty principle [ 13] , is 
readily obtained from a covariance phase diagram (Fig. 5). This is a plot of the covariance 
and the time derivative of the covariance as a function of time delay. The time delay 
corresponding to the first radial minimum of the phase diagram is defined as the 
coherence time. The rate of radial decay is a measure of the damping of the covariance 
and is directly related to the bandwidth-frequency ratio. That is, the phase diagram shows 
that the greater the bandwidth-frequency ratio, the shorter the piecelength required to 
obtain a wave packet description. We observe that the coherency time for a periodic 
process is infinite. In practice, some data sets must be truncated; therefore, the maximum 
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a. DATA COLLECTION SYSTEM b. TEMPORAL HISTORY OF 

INTENSITY FLUCTUATIONS 



COVARIANCE 



C. RESULTING CROSS-COVARIANCE CURVE 

Figure 4. A one-event, two-dimensional physical model. 
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piecelength will be set at four times the period associated with the mean frequency, 
unless the coherence time is less than half that value. In that case, the piecelength can be 
reduced to twice the coherency time without an information distortion. 
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D. Statistical Error and Period of Statioharity 

To facilitate the evaluation of the covariance curve and to establish the stability 
of the data set, it is necessary to utilize a set of statistical criteria. These criteria will now 
be developed. 

The sequential covariance curve is the basis' for testing the stability of the 
boundary conditions. The first class of stability test is the statistical error of the 
accumulative covariance that is produced as a result of the deviations of the covariance 
coefficient, at a particular time delay, between the different pieces of record. The 
statistical error is 
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where tp is the “student’s t’’ coefficient for a probability level, P, 

with m pieces [15]. The normal distribution would be used when there are more than 
30 pieces of data record. 

Briefly, let us examine the significance of this criterion with regard to the 
evaluation of the stability of the covariance curve. Assume that we have a set of 
piecewise covariance curves which contain some undesirable signal (noise). When these 
piecewise covariance curves are averaged together, we obtain the quantized accumulative 
covariance curve, together with the confidence limits, as shown in Figure 6. Only those 
parts of the covariance curve that are greater than the confidence limits are statistically 
significant to the probability level selected. Equations (32) and (33) are providing an 
estimate of the reproducibility of a covariance coefficient at a particular time delay. 
Since the noise has a low reproducibility compared to the information bearing signal, the 
noise lies below the confidence level and is thereby rejected. 

A second type of stability test is based on averaging the statistical error, from the 
accumulative averages, over the time delay range. The accumulative error. 
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QUANTIZED COVARIANCE 



implies that the average piecewise variance of the piecewise covariance curve is converging 
at the rate of 1/m with the variance of the accumulative covariance curve. This 
convergence can also be interpreted in terms of the bandwidth in the following manner. 
The variance about the mean of the piecewise covariance curve can be written in terms of 
the power spectrum as 
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where we are assuming a boxcar spectral distribution where all frequencies have the same 
amplitude. The limits of integration of this integral can be reduced to the effective 
bandwidth, . This means that there are just realizations per piece. Since the 

data record is T long, there are T realizations resulting from averaging. Thus, the rate of 
the convergence between the variance of the piecewise covariance and the accumulative 
covariance for stationary data is l/Bg^^T . In terms of the accumulative error, which is 

the square root of the variance, the rate of convergence is l/\/BgffT . When the 

accumulative error is plotted as a function of the inverse of the square root of integration 
time, we should obtain a straight line for a stationary process, whose slope is the square 
root of the inverse effective bandwidth (Fig. 7). This has been demonstrated rigorously 
by Jayroe and Su [11]. 

The accumulative error can also be plotted as a function of integration time to 
determine the integration time needed to retrieve the enhanced covariance curve. When 
the accumulative error stops decreasing as a function of integration time, then we are no 
longer enhancing the accumulative averages. At that point, we can terminate the 
integration. The period of enhancement can be clearly distinguished in the cold jet data 
shown in Figure 8, where we observe that the integration of additional data did not 
significantly enhance the results after 3 seconds of integration. The bottom diagram in 
Figure 8 clearly shows how stationary the boundary conditions were in this experiment. 

An important question still remains with regard to stationarity: “When does the 
data record become too nonstationary to correlate?” The statistical nonstationarity of the 
data does not imply it is meaningless. Consider, for example, the atmospheric wind 
velocity. Generally, over a 24-hour period the wind velocity is nonstationary, but for 
short intervals of time it is stationary. Thus, it is desirable to know when to create a new 
ensemble, since we are interested in isolating changes in environmental conditions that 
would result in a nonstationary data record. The chi-square test enables us to establish 
the limits of stationarity for the accumulative error [10]. The limits of stationarity are 
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The envelope defined by these chi-square limits (Fig. 9) sets the region of stationarity 
such that, when the accumulative error is within this envelope, the data record can be 
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ACCUMULATIVE ERROR <A<f» 



INVERSE OF THE SQUARE ROOT OF THE 

INTEGRATION TIME , sec"'''^ 

Figure 7. Accumulative error curve. 


considered analytically stationary. It is important to recognize that this is simply a 
restriction on the amount of variation of the bandwidth that can be tolerated. 
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Figure 9. Accumulative error curve for a nonstationary process. 


E. Probability of Retrieving Information About Transport Phenomena 

Utilizing the statistical tests that were just discussed, criteria will be established 
whereby an estimate can be obtained for the probability of retrieving significant transit 
time information concerning the local transport phenomena. 

The retrieval of transport information requires the detection of the common 
information contained in the temporal fluctuations from two spatially separated probes. 
The criterion of commonality requires the two signals to always maintain the same sign 
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over the period of integration. The probability of transport retrieval is thus identical with 
the probability that the accumulative covariance exhibits commonality and exceeds its 
statistical error. The degree of commonality between the average product can be treated 
as student’s t distributed if the data set is stationary, which is to say that this distribution 
apphes only to processes which belong to the same population as indicated by the 
chi-squared test. Under these conditions, the probability of retrieving the transport 
information can be computed directly from the student’s t distribution [16]. 

The degree of commonality is obtained for the student’s t from the ratio of the 
accumulative covariance and its statistical error, which is 
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The probability of commonality associated with this student’s t is then given by 
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dx , 


(38) 


where F j- is the gamma function [17]. The probability of commonality can now be 

plotted as a function of time delay (Fig. 10) and subdivided in accordance with geometric 
constraints into an information region and a noise region. The commonality probability 
because of noise, P^oise ’ maximum probability that exists in the noise region. The 

probability of transport retrieval, P^rans > 
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where Pjj^f is the maximum information commonality probability. The probability of 
transport retrieval must be 


P^ 2 50% (40) 

trans 
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PROBABILITY OF COMMONALITY 



THESE REGIONS ARE DEFINED BY THE EXPERIMENTAL CONSTRAINTS 


Figure 10. Probability of commonality. 


for the transport signature to be significant. This now affords an interval estimator for 
selection of the significant transit time signature by providing the probability that the 
common information from the two probes is obtained from the same streamline. 
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SECTION III. WAVE DESCRIPTION OF THE TRANSPORT PHENOMENA 


A. Introduction 

Experiments using correlation analysis [2-5,1 8] clearly show that the covariance 
curve contains signatures of the transport phenomenon, which cannot be retrieved or 
interpreted using the traditional time-series analysis. To extend the analysis, it is 
necessary to augment these analytical techniques with a new three-parameter wave 
description of the covariance curve. 

This description enables us to determine how the parameters of the wave packet 
are described in the covariance curve and thereby to develop a systematic classification 
procedure for the different kinds of turbulent transport processes. The wave description 
permits an interpretation of the covariance curve which separates the transport 
phenomenon into a convective component and dispersion component. 


B. Wave Description of the Covariance Curve 

The wave description will be obtained from a new wave packet algorithm, which 
will provide a three-parameter description of the signature associated with the 
fluctuations in the fluid transport. In the crossed-beam correlation technique, fluctuations 
are essential, in that their signature affords a time-dependent information element that 
enables us to retrieve the convection kinematics. The frequency description of the 
covariance delineations provides a field description of the affluence of the fluctuations 
present in the wave packet which, as was shown in Section I, are the signature of the 
dispersive kinematics. Before developing the wave description, the limitations of the 
Fourier frequency spectrum with regard to transport analysis will be considered. 

The covariance curve can be transformed from the time domain to the frequency 
domain using the Fourier analysis procedures outlined in the appendix. To demonstrate 
the effectiveness of this transformation, we can construct a covariance curve by the 
superposition theorem, which represents a wave packet, so that we know the input 
spectrum. For this example, suppose we obtain the covariance curve of our wave packet 
by superimposing 21 cosine functions of equal amplitudes which are equal frequency 
intervals apart over the range of 0.5 to 1.5 Hz (Fig. 1 1). The resulting covariance curve, 
shown in Figure 12, is damped because of the destructive and constructive interference 
introduced by superimposing the cosine functions. The energy density spectrum shown in 
Figure 13 is obtained when this covariance is transformed into the frequency domain. On 
inspection of the spectrum, we observe that, while the bandwidth of the output spectrum 
is identical with the input spectrum, the amplitude distribution of the output spectrum is 
not constant. This spectrum is correct. 
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Figure 1 1 . Graphical relation between input spectrum 
and energy density spectrum. 



Figure 12. Autocovariance curve in time domain. 


The explanation is as follows. The gain of the input energy density window was 
4.77 percent for each of the 21 frequencies. The entire bandwidth of these data occupies 
only nine wave numbers in our energy space, implying a mean gain of 11.11 percent. 
Figure 1 1 shows the input in terms of the output spectrum obtained. This discrete 
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PARAMETERS: 

AVERAGE FREQUENCY = 1.0 hz 

BANDWIDTH = 1.0 hz 

NUMBER OF WAVES -- 21 

FREQUENCY INTERVAL = 0.05 hz 



Figure 1 3. Energy density spectrum of a narrowband process. 
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spacing between the wave numbers, which is required for an orthonormal spectrum, 
causes an overlap and results in a distortion of the intensity distribution. (The wave 
number is that integer multiple of the inverse of the sampling interval that gives the 
frequency.) If there were integral multiples of wave numbers between frequencies, this 
would not happen. Thus, the price paid for obtaining the orthonormal frequency space is 
a loss in the resolution of the input spectrum. From this demonstration, we can conclude 
that the reproducibility of the actual frequency domain from the energy density 
spectrum, without either a foreknowledge of the desired results or a focusing of the 
spectrum, can be misleading. This forces us to seek a new parametric description that 
ehminates the masking effects of a Fourier analysis which resulted from the degree of 
uncertainty due to frequency quantization. 

This new frequency description of the time domain will be reduced to a 
description of the size, density, and location of an average envelope about a wave packet. 
Since the frequency domain is subject to distortion, our analysis is performed in the time 
domain. An algorithm for this wave packet description will be generated as a function of 
the average frequency <f> , the bandwidth B , and the number of frequencies N present 
(Fig. 14) in the wave packet. The superposition of N frequencies equally spaced and 
with the same amplitude, which is a boxcar distribution, is given by (for N-odd): 
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and 


Aco. = <a; > - u) . . 
1 1 


(44) 


AMPLITUDE 



Figure 14. Spectrum of the wave packet algorithm. 


Next, the cosine relationship for the difference and sum of angles gives 
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Simplifying the above relationship, we obtain 
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Now consider the summation in the above relationship: 
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then, for the (k + l)th term 
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That is, only the first and last terms in the summation remain, or 
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Returning to equation (47) and using the above results, we obtain 


(51) 
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(53) 



This algorithm provides a wave description of a packet of sinusoidal waves as a function 
of the average frequency, bandwidth, and number density. This can be interpreted, in 
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accordance with the superposition theorem, as an autocovariance curve that is a time 
domain description of the wave packet; that is, 





(55) 


Since this autocovariance curve is the result of the wave description in the frequency 
domain being transformed to the time domain, we shall refer to it as the wave packet 
algorithm. This algorithm can only be used to analyze a single wave packet description in 
the transport process. Where several different wave packets are occurring simultaneously, 
frequency discrimination techniques must be utilized to isolate the wave packet 
associated with each process. This relation is the wave description for the signatures of 
those fluid kinematics associated with the parametric fluctuations resulting from the 
dispersive transport process. The parameters (<f>,B, N) of this wave description are 
retrieved by the use of the least-squares fit of the wave packet algorithm to the 
autocovariance curve which is obtained from the suppression, of the transit time 
information in the cross-covariance delineation of two spatially separate probes. 

A more general form for the wave packet algorithm compresses the original three 
parameters into two, by redefining the time axis such that 


9 = <f> T 
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and defining the bandwidth ratio to be 
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These will be referred to as the similarity relations. Then, by this similarity transform, 
the wave packet algorithm, equation (55), becomes 
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which is the same as if the frequency were 1 Hz. 

A computer comparison of the wave packet algorithm with the equivalent 
description obtained from the superposition theorem shows that more precise results were 
obtained with the algorithm, especially when the number density of sinusoidal functions 
was high. The discrepancy between the two descriptions (maximum of 2.19 percent) 
which occurred for small covariance coefficients in the wings of the curve was a result of 
truncation errors induced by summing large numbers of sinusoidal functions. It can be 
concluded from this analysis that the wave packet algorithm is not only a faster but a 
more sensitive means of attaining the wave description than by using superposition 
techniques. 


C. The Wave Packet Interpretation for the Covariance Curve 

The parametric effects of the wave packet on the shape of the covariance curve 
will now be considered. 

Statistical processes are commonly grouped into three classifications: periodic, 
narrowband, and broadband processes [8] . Typical examples of the autocovariance curves 
associated with these processes, using data obtained from glow discharge experiments [4] , 
are shown in Figure 15. The distinguishing characteristic between these different processes 
• is the degree of damping associated with each process. 

To interpret the mechanism responsible for the damping, the wave packet 
algorithm will be reexamined. We observe that, when the algorithm is written in the form 


B, (27T<f>T) D(t) , 
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where 
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Figure 15. Autocorrelation curves for moving stria tions. 


we obtain a periodic cosine function of the average frequency and a time delay dependent 
factor, D(r) , which is clearly the damping factor in the wave description. 

Now consider the limiting case where the bandwidth of the wave packet is zero. 
Then the damping factor using THopital’s rule is 





( 61 ) 


for all time delays. This results in a periodic covariance curve of the average frequency, as 
we would expect. We can further conclude from an examination of the damping factor 
that the damping of the average frequency is a result of the number density and the 
bandwidth. In fact, to a first-order approximation, we can attribute the damping 
primarily to the bandwidth. 

Obviously the damping factor is a suppressing agent of the cosine function of the 
average frequency such that the bounds for the cosine function are set by the damping 
envelope. Using the same parameters in the wave packet algorithm that were employed 
with the superposition theorem (Fig. 11) in generating the covariance curve (Fig. 12), we 
can generate the autocovariance curve in Figure 16. This covariance curve generated with 
the algorithm has an important analysis impact in that it affords us the opportunity to 
examine the effects of the damping envelope. By inspection of this covariance curve, we 
can observe that the damping envelope defines the level of suppression for the cosine 
function of the average frequency. 

It is desirable to be able to obtain, by inspection, the average frequency of the 
wave packet directly from the auto covariance curve, as we would, for example, in the 
case of a periodic function. From an analysis of the wave packet algorithm, it follows 
that our freedom to retrieve the average frequency by inspection becomes more restricted 
as the bandwidth increases, because it will result in an increase in the suppression of the 
characteristics of the average periodicity. In fact, there is a point where the bandwidth 
suppression of the covariance is so strong that even the first quarter cycle of the average 
periodicity is suppressed, thereby preventing the retrieval by inspection of the average 
frequency associated with the wave packet. 

It is thus necessary to establish a classification system which divides the statistical 
processes into those processes that are amenable to a direct information retrieval scheme 
and those processes that require a more sophisticated information retrieval scheme. This 
is the objective of the next section. 


D. Classification of Covariance Curves 

A procedure for the systematic classification of covariance delineations is 
developed that enables us to group statistical processes in accordance with the minimum 
data processing required for the retrieval of the wave packet description. 
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There is a classical system defining the classification of statistical processes, as 
pointed out in the last section. No damping implies a periodic process, some damping 
imphes a narrowband process, and strong damping implies a broadband process — thus a 
heuristic definition. While this definition is an aid in grouping data, it lacks an analytic 
formalism. More specific definitions of the processes are to be found in optics [13] , but 
they do not directly apply to our discussion. 

The behavior of the damping envelope, as pointed out in the last section, is our 
primary concern. The principal factor governing the damping envelope is the bandwidth 
of the statistical process. From the similarity transform of the wave packet algorithm 
[equation (58)] , we observe that a behavior pattern for the damping envelope also can be 
obtained in terms of the bandwidth ratio [equation (57)] . The beauty of the 
dimensionless bandwidth ratio lies in its ability to render a generalized damping model. 
Since we would like our classification system to have significance in the frequency 
domain as well as the time domain, we will initially direct our attention to the behavior 
of the spectral envelope as a function of the bandwidth ratio. 

The range of the bandwidth ratio is from zero to two for the defined boxcar 
spectral envelope. If the bandwidth exceeds two, then the spectrum folds (Fig. 17) 
because of the even property associated with cosine functions. It follows that the three 
analytical classifications, from the spectral standpoint, of the correlation curve based on 
the bandwidth ratio, ^ , are: 

1. Periodic processes: 0 = 0 . 

2. Narrowband processes: 0 < 0 < 2 . 

3. Broadband processes: 0 > 2 . 

With this proposed classification system as a guide, a closer examination of the time and 
frequency domains is in order. 

To facilitate this examination, the computer program, whose flow diagram is given 
in Figure 18, was used to generate a set of covariance curves to demonstrate the effects 
of the bandwidth ratio on the damping envelope. Our objective is to illustrate that the 
classification system just defined not only has significance in the frequency domain, but 
that this system also groups covariance curves such that in periodic and narrowband 
processes the average frequency can be retrieved by inspection in accordance with 


<f> 


4 T 


FCP 


(62) 


37 



A. NARROWBAND LIMIT (j8 = 2) 


GAIN 



B. BROADBAND ( > 2 ) 


GAIN 



Figure 17. Schematic representation of the envelope for the 
spectrum from narrowband and broadband processes. 


where is the delay time of the first crossing point which is taken as the delay 

corresponding to the initial zero value of the covariance coefficient. This relation 
[equation (62)] indicates that the first quarter cycle of the covariance curve corresponds 
with the first quarter cycle of a cosine wave of the average frequency. 
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Figure 18. Flow diagram of spectral decomposition program 








The results for mappings of the autocovariance curve over the bandwidth ratio 
range of zero to eight, of which three typical results are given in Figure 19, clearly 
indicate that the average frequency can be retrieved by inspection in periodic and 
narrowband processes. However, in broadband processes, where the bandwidth ratio 
exceeds two, the damping envelope always suppresses the first quarter cycle of the 
average frequency, making information inspection retrieval impossible. This behavior is in 
accordance with our original classification model objectives in which we wanted to 
differentiate between the statistical processes where the average frequency could be 
retrieved by inspection from those processes where analytical techniques are required. 

Analytical techniques for ascertaining the wave packet parameters include spectral 
decomposition (Appendix) and regression analysis of the wave packet algorithm. Link 2 
of the computer program (Fig. 18) was used to test the information retrieval capabilities 
of the energy density spectrum. The covariance curves subjected to this spectral 
decomposition were generated by the wave packet algorithm. This means that the output 
spectrum should be identical to the boxcar input spectrum (Fig. 14). The boxcar 
spectrum is obtained when there is an integral wave number spacing between each of the 
input spectral frequencies (Fig. 20). This is a special case, and in general we would not 
have this integral spacing. Therefore, it is much more difficult to obtain the wave packet 
description directly from the spectrum. To demonstrate this point, observe the 
narrowband spectrum shown in Figure 21. In this example, the quantization effect caused 
by discrete sampling has suppressed the information regarding the distribution and the 
number density. It can be shown that the description of the time domain in the 
frequency domain is limited to the same number of unique representations, as pointed 
out in the explanation of the spectrum in Figure 12. Thus, in addition to the distortion 
illustrated in this figure, there is also an uncertainty associated with the spectrum which 
is equal to the wave number spacing. This uncertainty in the spectral description renders 
the spectrum unsuitable for the detection of small changes in the average frequencies 
between different covariance curves. Therefore, the application of change statistics of the 
average frequency for periodic and narrowband processes requires either inspection 
retrieval or wave packet regression analysis. 

In broadband processes, the Fourier spectral description of the autocovariance 
curve will deteriorate. Using the same spectral input model as used in the narrowband 
discussion, but with a broader bandwidth, we obtain the integral wave number spectrum 
shown in Figure 22. We would expect a result similar to that shown in Figure 20, except 
it would be folded. Instead, we get a spectrum with a marked ringing signature in the 
folded region. The deterioration results from the fact that the Fourier spectrum, which is 
a correlation between the covariance curve and a sinusoidal function, is being calculated 
for a covariance curve which no longer has a dominant signature of the average 
frequency. This signature has been suppressed by the bandwidth damping. Naturally, the 
broadband nonintegral wave number spectrum of the auto covariance curve (Fig. 19c) 
results in even more deterioration for the output spectrum (Fig. 23). Unlike the 
nonintegral narrowband distribution, the average frequency of the wave packet will not 
tend to dominate the broadband spectrum, and must be retrieved in terms of a model. In 
fact, our whole analysis of a broadband process requires a model to interpret the wave 
packet parameters. 
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Figure 19. Autocovariance curves of narrowband 
and broadband processes. 
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Figure 20. Transform of square input with integral spacing. 
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To verify the validity of the wave packet description, we must examine both the 
phase diagram (Fig. 5) and the accumulative error (Fig. 7). The phase diagram affords an 
independent test for the average frequency of the wave packet, since the derivative is 
normalized by the average frequency. The effect of this normalization is to remove the 
oscillation in the radial vector resulting from the average frequency; therefore, only the 
proper frequency normalization removes the oscillations. The average frequency defined 
by the wave packet description is compatible with the restrictions imposed by the phase 
diagram. The effective bandwidth (Fig. 7) obtained from the accumulative error can be 
shown to be the equivalent of the bandwidth in the wave packet description. It can 
therefore be concluded that the wave packet description is compatible with our overall 
statistical model. Thus, it is a vaUd model for interpreting the spectra from any process. 

Since our model for interpreting the spectrum is vahd, our problem in broadband 
analysis is again the quantization uncertainty introduced by discrete data. To circumvent 
this problem, which is inherent to the frequency spectrum, the least-squares criterion for 
fitting the wave packet algorithm to the autocovariance curve was used. To demonstrate 
the effectiveness of this criterion, the residue error mapping of the algorithm fit to the 
covariance curve was made. A typical example of the residue mapping from fitting to the 
algorithm to covariance curves of boxcar spectral inputs is shown in Table 1. The wave 
packet algorithm retrieved the exact parameters for the average frequency, bandwidth, 
and number density in all cases. When the narrowband process was simulated using an 
elliptical input spectral window, the exact average frequency was retrieved and the 
bandwidth retrieved gave an equivalent area under the spectral curve (Fig. 24). Similar 
results were obtained for broadband processes (Fig. 25) using the algorithm. The greatest 
error in the parameters of the wave number algorithm was found in the number density, 
where a unique value was not always clearly determined. In spite of this limitation of the 
wave packet description, the results were always reproducible with the original covariance 
curve to less than 2.61-percent deviation for narrowband processes and less than 
3.46-percent deviation for broadband processes. In general, the deviation between the 
known covariance curve and the one generated by the algorithm parameters occurred in 
the far wings of the curve. The number density can be viewed as a fine-tuning 
adjustment, since it only shghtly affects the damping function. A comparison of results 
obtained with the wave packet algorithm provides the sensitivity and dependability 
required for transport analysis. 

Thus, the classification system which has been just developed not only is an 
analytical guide to the statistical techniques apphcable to the analysis, but it also affords 
an index to the spectral type. While Fourier analysis affords a parametric estimate for the 
wave packet description, it lacks the necessary precision for change detection which is 
obtained with either the direct inspection retrieval or the wave packet algorithm. A 
degree of care must be exercised with the direct retrieval of the average frequency 
because its apphcation is Umited to periodic and narrowband processes. 

With this basis for the retrieval of the wave packet information, we can now 
examine the analysis of the transport process. 
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TABLE 1. MAPPING OF ALGORITHM DEVIATION RESIDUE 
FOR AN AVERAGE FREQUENCY OF 5 Hz 


NUMBER 
OF WAVES 

BANDWIDTH (Hz) 

5.400 

5. 600 

5.800 

6.000 

6.200 

6.400 

6.600 

6. 800 

22 

0.00312 

0.00217 

0. 00091 

0. 00050 

0.00174 

0.00273 

0.00331 

0.00363 

24 

0.00316 

0.00225 

0.00103 

0. 00033 

0.00158 

0.00259 

0.00321 

0.00356 

26 

0.00320 

0.00232 

0.00113 

0.00020 

0. 00144 

0.00247 

0.00313 

0.00351 

28 

0.00323 

0.00237 

0.00121 

0.00009 

0.00133 

0.00237 

0.00306 

0.00346 

30 

0.00326 

0.00242 

0.00129 

0.00000 

0.00124 

0.00229 

0.00300 

0. 00342 

32 

0.00328 

0.00247 

0.00135 

0.00008 

0.00116 

0.00222 

0.00295 

0.00339 

34 

0. 00331 

0.00250 

0. 00141 

0.00015 

0.00109 

0.00215 

0. 00290 

0.00335 

36 

0.00333 

0.00254 

0.00146 

0.00021 

0.00103 

0.00219 

,0.00286 

0.00333 

38 

0.00335 

0.00257 

0.00150 

0.00026 

0.00097 

0.00205 

0.00283 

0.00330 

40 

0.00337 

0.00259 

0.00154 

0. 00031 

0.00093 

0.00200 

0.00279 

0.00328 

42 

0.00338 

0.00262 

0.00158 

0.00036 

0.00088 

0.00196 

0.00276 

0.00326 




STANDARD DEVIATION BETWEEN ALGORITHM 






AND TEST CURVE 





<f> = 5 Hz 

Parameters of Test Curve: B = 6.0 Hz 

N =30 waves 


E. Separation of Convective and Dispersive Transport Phenomena 

The various statistical techniques that have been discussed will now be unified to 
obtain a method of separating the convective and dispersive transport information from 
the signatures contained in the covariance curve. 

Assume that the temporal histories of the fluctuations in the optical density of a 
flowing fluid have been retrieved at three points along the streamline (Fig. 26). Three 
cross-covariance curves can be obtained by cross-correlating the temporal information 
obtained from the x-probe with the temporal information retrieved from the three 
yj-probes. ' This set of three curves provides a mapping along the streamhne of the 

time-averaged transport. To establish the validity of this transport description, the 
accumulative error is examined to determine whether the process is stationary. If we have 
a stationary ensemble, then peak identification procedures are used to determine the 
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Figure 24. Wave packet algorithm fit of elliptical envelope. 


probability of retrieving the transport and of isolating the transit time of the transport. 
The two transit times thus retrieved are plotted as a function of their transit distance. 
The slope of the straight line passing through these points is the average convection 
velocity along the streamline. 

To obtain the dispersive signatures, it is necessary to transform the covariance 
curve from a fixed representation to a representation moving at the convective velocity. 
Two techniques can be used to obtain this transformation. The first is to suppress the 
convective information by transforming the gain spectrum into a fourth-order covariance 
curve. The second technique is to directly obtain a fourth-order covariance curve of the 
cross-covariance curve by the relation 


(xy)2= 

T m 



T -Tm T +Tm 

_£ e 

2 


(63) 
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GAIN 



where is the transit time. This second suppression technique will be used in these 

discussions because of its simplicity. (The computer program in Figure 18 is designed to 
use the first technique.) 

This fourth-order covariance curve contains only dispersive information since this 
suppression technique has effectively suppressed the transit time information concerning 
the convective process. The dispersive information is retrieved from the fourth-order 
covariance by fitting the wave packet algorithm to it and obtaining the parameters for 
the wave description. 
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Figure 26. Retrieval of convective and dispersive transport. 
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The dispersive velocity, v^j , can be expressed in terms of the change in the 
average frequency of the wave packet, the convective velocity, v^ , and the transit 

distance. The following argument is used. The change in the packet size, As, due to 
dispersion is 


As = - 






(64) 


where the quarter cycle spread represents the change in the packet shape. The time 
interval, , over which this change occurs is 


T 

P 



(65) 


It follows then from the elementary laws of kinematics that the dispersive velocity is 


^ 4 <fapxfap> 


( 66 ) 


The convective and dispersive velocities can thus be separated by first retrieving 
the convection from the second-order cross-covariance curve, then suppressing the 
convective information and determining the wave description of the fourth-order 
autocovariance. The change of average frequency as a function of transit distance affords 
the dispersion. Returning to the continuity equation [equation (9)], the convection is the 
mean transport; whereas, the dispersion is the turbulent flux. 
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SECTION IV. APPLICATIONS TO REMOTELY SENSED DATA 


A. Introduction 

The empirical evidence of applicability of the statistical techniques, which have 
been discussed, to remotely retrieve the transport kinematics will now be demonstrated. 
Since heuristic arguments regarding the crossed-beam test arrangement have thus far been 
presented, we begin by presenting the analytical arguments that demonstrate the 
crossed-beam correlation technique’s potential to remotely utilize detected optical 
modulations to obtain local information regarding the turbulent transport process. 

The applicabihty of the kinematic description obtained with these techniques will 
be illustrated under three different environmental conditions. The analysis of moving 
striations in the glow discharge is an example of one-dimensional flow which illustrates 
ensemble grouping, process identification, and transport analysis. The analysis of the 
convective and dispersive transport associated with the plume of an air jet is used to 
show the applicability of these techniques to a stationary, two-dimensional turbulent flow 
in a fixed direction. The third illustration involves transport phenomena in the 
troposphere where the turbulence flow is nonstationary, and, because the flow direction 
is a variable, it is also three dimensional. 


B. The Crossed-Beam Test Arrangement 

The fundamental crossed-beam correlation method which was initially conceived 
by Krause et al. [1] can be perceived most readily by examining Figure 27. Assume that 
a region of turbulent flow is contained within the broken line and is being convected in 
the direction normal to the plane of the diagram. Two optical systems consisting of a 
radiation source and a radiometer are arranged so that collimated beams of radiation pass 
across the flow in two mutually perpendicular directions and the beams intersect at the 
point to be investigated. The wavelength of this radiation is selected so that it is 
modulated by one or more species of the flow, allowing fluctuations of the concentration 
or density of the selected species to be reflected in changes in the intensity of the 
radiometer [ 19] . Each beam alone reflects only an integral of the appropriate 
fluctuations along its entire path length; however, it will be shown that the covariance of 
the signals from the radiometers does yield local information about the point of 
investigation. 

The retrieval of this local information can be explained intuitively as follows; The 
instantaneous signal at each detector is the resultant of all modulating influences along 
the path of the beam at a particular time. These modulations of the source radiation can 
be categorized into two groups: 


52 




Figure 27. Schematic cross section for the crossed-beam 


monitoring system in a turbulent flow. 


1. Those which occur sufficiently close to the point of beam intersection to 
introduce a related (i.e., correlated) modulation in both beams. 

2. The remaining modulations that occur at a sufficient distance from this 
point and are uncorrelated because they introduce an unrelated effect on the beam 
intensities. 

The correlation of the two signals will yield an average value of zero for those 
portions of the modulated signal created by the unrelated flow fluctuations; whereas, the 
related flow fluctuations will yield a finite average product. In principle, then, the 
measured covariance is a function of only those fluctuations that modulate the beams in 
the correlation area around the beam’s intersection point. 

To analytically demonstrate that the covariance of the two detected modulations 
yields a measure of required turbulent properties, we will introduce the following 
coordinate system. The point of beam intersection has the coordinates (x, y, z) where 
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the X and y axes are oriented parallel to beams one, Si Ri , and two, S 2 Rj , respectively. 
The distances measured from the point of beam intersection are denoted 
by f , 7? , and | in the x , y , and z directions, respectively. 

For beam one, the intensity of the radiation recorded at detector Ri at 
time t is 




lo exp 



K(x + y. 


z, t) df 


(67) 


where lo denotes the initial intensity of the beam and K is the appropriate extinction 
coefficient. Here, “extinction" is any possible mode of beam attentuation. Two common 
forms of extinction are pure absorption by a flow constituent and scattering by 
particulate matter in the flow. The extinction coefficient, K , is also a function of both 
the wavelength of the radiation and of the thermodynamic parameters of the 
constituents [ 19] . Since this dependence does not add enlightenment to the present 
discussion, it will not be exphcitly shown. 

Whatever the actual mechanism of extinction, it should be selected such that the 
extinction depends on a required flow property. Since the flow properties are functions 
of both position and time in a turbulent flow, the extinction coefficient will be similarly 
dependent, and thus we can just refer to fluctuations or modulations of the extinction 
coefficient. Since these changes are linearly related to fluctuations of a flow property, 
statistical properties of the flow do represent statistical properties of the extinction 
coefficient. 

The instantaneous extinction coefficient, K , can be written as 


K(x + y, z, t) = k(x + y, z, t) + K(x + y, z) 


( 68 ) 


where k is the instantaneous fluctuation of the extinction coefficient about the mean 
value, K . The detected intensity [equation (68)] can now be written as 


Ii(t) = exp 


Ri ) ( 

/ K(x+ y, z) > exp< 
Si . ) ( 


- / k(x + y. z. t) 
Si 


(69) 
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If the extinction process is such that the integral of the fluctuations [/k(x + 
f , y, z, t)df] is sufficiently small to permit linearization of that portion of the exponential, 
using the series expansion of an exponential to first-order approximation, the intensities 
can be written as 


Ii(t) 


Ifl exp 



K(x + y, z) 



X 


Rl 

1 - / k(x - 1 - y, z, t) 

Si 


(70) 


This assumption does not restrict the crossed-beam correlation technique to applications 
with small modulation of the optical beam because of the following: 

1. The integral in question represents a sum of a number of statistically 
independent events, which will tend to reduce the value of the integral. 

2. It was shown [20] that if the integral of the fluctuations is of the order of 
10 percent of the mean integrated value or less, then an optimum value for the mean 
attenuation is given by 

Rl 

f K(x + £, y, z) df = 1 . (71) 

Si 


For this magnitude of fluctuation, the linearization would be acceptably accurate. 

3. If larger fluctuations relative to the mean value are experienced, it would be 
acceptable and indeed desirable to reduce the mean absorption, by reducing the gain, so 
that linearization is still possible. 

If the intensity at the detector is now written as the sum of its time 
average Ii and a fluctuation component i, (t) relative to the average value, it can be 
shown that 


_ Rl 

ii(t) = - ^1 J k(x + y, z, t) d£ . (72) 

Si 
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Within the limits of the above discussion, we obtain the expected result that the 
fluctuation of the intensity at the detector is proportional to the instantaneous integral 
of the modulations of intensity along the entire light path. 

By inspection, we can write similar results for the other beam: 


_ ^2 

i2(t) = - I2 / k(x, y + 7 ], z, t) d?} . ( 73 ) 

^2 


If we now take these two fluctuating intensities, i,(t) and i^Ct), and determine their 
time-average product, we obtain the covariance in a piecewise form similar to 
equation ( 24 ); i.e., 

^ J J T / k(x + ^, y, z, t)k(x, y + T), z, t) dt dn , 

Si S 2 0 (74) 


where the spatial and temporal integration has been reversed. This now gives us the 
covariance as a function of the modulations of the extinction coefficients, or in terms of 
the statistical flow parameters of the fluid. 

To most conveniently understand how spatial resolution is obtained with the 
cross-covariance, we will initially consider just the temporal integration at 7 = 0: 


p(x + y + T], z) 


= i f k(x+ 4 , y, z, t)k(x, y+T], z, t) dt. 

T=0 0 

( 75 ) 


This term clearly represents the covariance of the fluctuations occurring at the two points 
(x + f,y, z, t) and (x, y + 17, z, t). If one or both of these points are sufficiently far 
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from the beam intersection point in an inhomogeneous flow, the fluctuations will be 
mutually random and the resulting covariance will be zero. Only those points contained 
within the correlation area around the beam intersection point will contribute to the 

measured value of the covariance. Therefore, the value of the covariance, does not 

change if these limits are replaced by those limits of the locally correlated area, although 
formally the limits of spatial integration in equation (74) extend from source to detector. 
Because of our prerogative in establishing the limits, we can view the measured quantity 
as reflecting only the local modulating information and thus providing spatial resolution. 

To retrieve the transit time between beams, the two beams must be separated in 
the z-direction a distance | (Fig. 28). To the flow, there still appears to be an 
intersection point for the beams, which we shall refer to as a “virtual crossing point.” 
TTiis virtual point is the only point at which a particle moving in the z-direction can pass 
through both beams; therefore, the intersecting beam argument for spatial resolution still 
holds except that the information is correlated now at a time delay. The time delay , 
corresponding to the maximum covariance, is a measure of the transit time between the 
beams. 



Figure 28. Crossed-beam test configuration for the cold jet. 
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This now provides an analytic association of the extinction coefficient with the 
crossed-beam correlation model. This model could be easily expanded to include 
spectroscopic emission and absorption [19], but for this discussion, our attention can be 
focused on just the kinematics of the flow. 


C. Ion and Electron Transport in the Glow Discharge 

To empirically demonstrate the feasibility of the retrieval and separation of 
convective and dispersive fluid transport, the application of the crossed-beam correlation 
technique to an investigation of moving striations in the flow discharge will be discussed. 

The striation is the disc of ionized and excited gas found in the positive column 
of the glow discharge (Fig. 29). Since the ions normally give the striated gas a net 
positive charge, these striations, which we will call forward-moving striations, move 
toward the cathode (Fig. 30) [21,22]. This emission process is the signature for the ion 
transport. 



Figure 29. Experimental setup for glow discharge. 
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There exists in the glow discharge a second transport; namely, an electron flow 
from the cathode to the anode. The electron flow is an extinction process whose 
signature modulates the ionization and excitation of the striation resulting in detectable 
undulations in the emission. When the striations are weakly ionized compared to the 
electron flow, the backward-moving striation exists (Fig. 31). Other investigators have also 
made this observation [23,24] . The importance of the backward-moving striation hes in 
the fact that it demonstrates the dispersive and convective transport characteristics of the 
electron flow. The glow discharge consequently provides an excellent stationary 
one-dimensional data source to demonstrate the effectiveness of the statistical method 
and algorithms used for fluid transport analysis. 

The transport process in the glow discharge was monitored by the crossed-beam 
configuration shown schematically in Figure 29. While a wide range of different kinds of 
transport was observed [4], our primary interest in these discussions is the data obtained 
when both electron and ion transport were simultaneously detected. A typical segment of 
such a data record is shown in Figure 32. 

The initial step in the data analysis was to cross-correlate the data and obtain the 
accumulative covariance, the accumulative error, and the transport probability in 
accordance with the data reduction scheme shown in Figure 33 [25]. The accumulative 
error was used to determine when the discharge was in equilibrium. It was found that 
about a half-hour warmup time was required for the transport process to reach its first 
stage of equihbrium. Once this stage was reached, the data could be treated as stationary 
until there was a discrete transition to the next lower mode. These transitions continued 
until the ground mode was reached. The physical significance of these modes will be 
explained. For now it is important to recognize that the accumulative error did permit us 
to differentiate between the different modes and to group coherently the data in its 
proper mode. 

The transport probability showed: 

1. A high probabihty that two different transport processes were occurring 
simultaneously in different directions — one a high frequency transport, the other a low 
frequency transport. Retroactive filtering [26] ‘ was used to separate these two transport 
processes. 

2. A high probability that the fluctuations in the discharge current were related 
to the transport moving toward the anode. Spatial mappings of the backward-moving 
striations showed there was indeed a correlation between these striations and the current 


1. Essentially, retroactive filtering is where the covariance curve is transformed to the 
frequency domain. The gain of the frequencies to be filtered are suppressed by setting 
them to zero and transforming the remaining energy density spectrum back to the time 
domain. Thus, the covariance curve is filtered without reprocessing the data. 
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Figure 32. Separation of electron and ion transport. 


fluctuations. Thus, the transport probability in these experiments served to isolate the 
different processes occurring in the glow discharge during a time interval which has been 
shown, using the accumulative error, to be stationary. 

Based on the transport probability, if we retroactively suppress the high frequency 
information by filtering the processed results, we obtain a field description of striations 
moving toward the cathode. The results of the forward-moving striations are summarized 
in Figure 34. The spatial mapping of the covariance curves affords a representation of the 
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Figure 34. Ion striation model. 


behavior of the ion striation, in terms of the average time, as a function of the transit 
distance between the optical beams. The transit time corresponding to the maximum 
significant covariance, when plotted as a function of the transit distance for a number of 
beam separations, can be fitted with a straight hne whose slope is the average convection 
velocity for the forward-moving striations. The positive slope of this line means, in 
accordance with the test configuration, that the transport is from the anode to the 
cathode. The convection velocity obtained for these striations ranged from 0.5 to 
10 km/sec. 
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The fourth-order autocovariance mapping in which the observer is in a moving 
frame where the transit time information has been suppressed indicates that the 
longitudinal length scale of the forward-moving striation was invariant as it was 
propagated along the positive column of the discharge. The proof of this invariance is 
obtained from the spatial independence of the wave description. The longitudinal length 
of the ion striation is obtained from the relation 


i 

X 


v 


c_ 

2<f> 


(76) 


where v^ is the convection velocity and <f> is the average frequency obtained from the 

wave packet algorithm. The experimental results for this striation length varied from 5.3 
to 10.2 cm as a function of the electrical energy supplied to the discharge. These length 
scales are in keeping with values obtained from high speed movies of the discharge. 

The repetition frequency, which is the rate at which an event occurs at a point, 
can be retrieved either directly from the crosscovariance or from the frequency spectrum 
of this curve by retroactively suppressing the wave packet description. The remaining 
energy density spectrum is due to the rapidity of the process. (If no spectrum exists, 
then the average frequency is the repetition frequency.) Using the repetition 
frequency, fj. , and the convection velocity, v^ , the longitudinal wavelength is 

V 

X = -r • (77) 

X f 

r 


The longitudinal wavelength of the ion striation was always an integral multiple of the 
length of the positive column, ; i.e., 

{. = n\ where n = 1, 2, 3, 4 . (78) 

pc X 

X 


The highest mode measured was four. The most stable mode was n = 1 and will be 
referred to as the ground mode. Experimental observations revealed that there was always 
a discrete transition between the modes, which could be detected by changes in the 
accumulative error. These changes appeared when the power supplied to the discharge 
was changed, whereupon a nonground mode would occur. This higher mode would decay 
back to the ground mode in discrete transitions. For standing striations, the wavelength 
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relation changed such that the ground mode corresponded to a wavelength equal to the 
distance between like points of adjacent striations rather than the length of the positive 
column. 

The signature of the electron transport was enhanced by suppressing the lower 
frequencies (Fig. 35). The analysis used for the ion striations was' repeated for the 
electron transport. Since the test configuration has been reversed, the positive convection 
velocity (Fig. 36), which is obtained from the plot of the transit time as a function of 
transit distance, implies that the flow direction is from the cathode to the anode. This 
flow was identified as an electron transport by first obtaining a cross-covariance between 
the optical fluctuations in the ionized gas and the fluctuations in the discharge current, 
and then calculating the transport probability. Close analysis of the temporal histories of 
these fluctuations revealed that there was not a continuous flow but rather a discrete 
transport; therefore, the electron transport can be more accurately referred to as a flow 
of electron packets which are propagated toward the anode. 

When the transit time information was suppressed [equation (63)] , the mapping 
of the fourth-order autocovariance curves (Fig. 36) showed that the electron packet’s 
shape was not invariant during its propagation toward the anode. This observation was 
supported by a change in the average frequency obtained from the wave packet 
algorithm. The electron packet’s dispersion velocity [equation (66)] was calculated to be 
approximately 260 m/sec compared to the ion convection velocity of 16 km/sec. Using 
the relation for the longitudinal length scale [equation (76)] , the electron packet size 
changed, on an average, from 1.0 to 1.5 cm over a transit distance of 30 cm. The major 
problem encountered in this work was that the dispersive effect required a resolution of 
one part in 10’ per centimeter t^o retrieve the changes in the average frequency. Thus, 
only by the use of averaging techniques and the wave packet algorithm was it possible to 
obtain the consistent results required for this change detection. 

These typical results demonstrate the feasibility of the analysis procedure used. 
Systematic analysis afforded by statistical techniques enhances the results and thereby 
reduces the amount of interpretation required. Since moving striations and turbulence 
eddies are both treated as moving packets of molecules, we can, in our next discussion of 
turbulent transport, conceptually approximate the turbulence eddies as being moving 
striations. 


D. Convective and Dispersive Transport in a Supersonic Jet Shear Layer 

The supersonic cold jet transport will be used as an empirical illustration of the 
applicability of the convective and dispersive transport model to stationary 
two-dimensional transport problems. To obtain a two-dimensional coordinate system, the 
axial streamline is assumed to be an axis of symmetry. This two-dimensional model is a 
direct extension of the one-dimensional description used with the glow discharge except 
that it accounts for lateral kinematics that could be neglected in the one-dimensional 
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Figure 35. Electron packet flow in the glow discharge. 
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flow using the constraints of the system. Our conceptual view must become more 
sophisticated for transport processes which are not subjected to the constraints of a 
one-dimensional flow. The shadowgraph of the cold jet flow clearly shows the flow 
expansion, which must be accounted for in any complete description of the kinematics. 
To this end, Corrsin’s spherical description is introduced to account for this transport 
behavior [27] . The applicabihty of the spherical model was demonstrated in the 
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discussion of the fluid model. The implications of this assumption are a discrete flow 
formed from a collection of expanding spheres emanating from a point source. These 
spheres are referred to as “turbulence eddies,” and are assumed to be of constant mass 
with a changing volume due to thermodynamic and force gradients. The current state of 
the statistical information retrieved prohibits a detailed consideration of these gradients, 
thus limiting the predictability of any of the algorithms to a local description. This means 
that the long-range kinematics can be extrapolated only from a spatial mapping. 

The cold jet is schematically shown in Figure 28. Again, we will assume a fan 
beam test arrangement, which is obtained by making a series of measurements along the 
streamline under similar test conditions. Using those data obtained along a streamline, by 
Fisher and Johnston [2] , the spatial mapping of the covariance curves shown in Figure 37 
was obtained. When Fisher and Johnston plotted the transit time as a function of transit 
distance, they obtained a convection velocity of 347 m/sec. They further showed that the 


TRANSIT DISTANCE , cm 
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data were stationary by using the accumulative error (Fig. 9) to establish the reliability of 
their results. The decay in the amplitude of the peaks of the covariance curves as a 
function of transit distance is a signature indicating the presence of a dispersive transport. 
To emphasize the need for a more sensitive wave description, we show the typical 
frequency spectrum of these covariance curves in Figure 38. Upon inspection of this 
spectrum, it is obvious that a resolution of one part in 10’ per centimeter for the change 
in the average frequency, which is required to retrieve the dispersive transport, cannot be 
obtained. Thus, we are again, as with the glow discharge, forced to use regression analysis 
of the wave packet algorithm [equation (55)] for each of these covariance curves to 
retrieve the parameters for the wave description. 

The bandwidth ratio did change as we moved radially outward across the shear 
layer from 1.47 to 1.92. While the plume of the jet involves a clearly narrowband 
process, we observe that, as we move out of the jet’s flow, the process tends to approach 
a broadband process. This is the expected result since atmospheric turbulence is 
characteristically a broadband process. 

Using the average frequency obtained with the wave packet algorithm, in the 
dispersion relation [equation (64)] , an average radial dispersion velocity of 19.7 m/sec 
was obtained with a convection velocity of 347 m/sec. If we assume that there is 
spherical expansion from a point source, in accordance with Corrsin’s model and our 
fluid model, then we can obtain a lateral mixing rate, L , given by 


L 



(79) 


Applying this relationship to this data set, we obtain a lateral mixing rate of 0.1 135. This 
rate implies a radial plume expansion of 11.35 cm/m downstream. This expansion rate is 
plotted on the shadowgraph of the jet plume, along with the locations of the points from 
where the local turbulent information was retrieved (Fig. 39). The outer plume cone is 
drawn so as to include all of the plume; whereas, the inner cylinder represents no plume 
expansion. It is clear from studying this shadowgraph that the cone defined by the rate 
of radial expansion is the average rate of expansion. This follows logically, as would be 
expected, since the source of the statistical information is obtained from the covariance 
curve, which is an averaging process, and is in accordance with the average turbulence 
flux [equation (12)] . Thus, the rate of dispersive transport obtained from the wave 
description affords an accurate analytical description of the physical situation associated 
with the cold jet plume. 

It is recognized that this prediction of convective and dispersive transport is only 
short range, since long range prediction would require an accounting for shear and 
gravitational forces [6,28,29] . In spite of this limitation, this new description of 
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Figure 39. Prediction of lateral scales in cold jet. 




turbulent transport is a significant contribution to turbulent investigations, because, for 
the first time, it accounts for the turbulent flux in the jet flow. 


E. Convection and Dispersion in the Atmosphere 

The crossed-beam correlation technique has been applied to the measurements of 
both clear air turbulence [3] and smoke stack emission.^ The experiments involving the 
retrieval of clear-air turbulence convection velocities and scales served to establish the 
statistical models for recognizing and interpreting nonstationary data. The geometric 
constraint imposed by using just two single-beam radiometers (Fig. 40) prohibited the 
application of change statistics which are required for the retrieval of the dispersive 
transport [30,31] . This hmitation was overcome in the smoke plume experiments by 
introducing a fan-beam radiometer which enables us to retrieve the local data at five 
different positions along the streamline. 


THE LIGHT SOURCE IS SCATTERED SOLAR RADIATION (PASSIVE SOURCE) 



THE FLUCTUATIONS IN LIGHT 
INTENSITY ARE RECORDED 


Figure 40. Atmospheric crossed-beam detection 
system for clear-air turbulence. 


2. Bilbro, J.W., et al.: Passive Remote Detection of Smoke Plume Height and Velocity 
using Correlation Techniques. To be published as a NASA Technical Memorandum. 
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To apply the statistical techniques discussed for cold jet experiments, we must 
first obtain a covariance curve which affords a stationary description of the transport 
process. In the jet plume and the glow discharge experiments, the test for stationarity 
was done primarily as an aesthetic step to obtain completeness. Such is not the case for 
atmospheric transport phenomena, because the interval of stationarity in the atmosphere 
is usually less than an hour in duration [32] . Therefore, the accumulative error must 
always be checked to determine whether it remains within the error bounds set by the 
chi-square test. Normally, the accumulative covariance is checked for stationarity and is 
used when it meets this requirement. If the accumulative error curve indicates that 
accumulative covariance is nonstationary, then we normalize each piecewise covariance 
curve and average them to obtain the quantized covariance. The accumulative error of 
this quantized covariance curve is then checked for stationarity to determine where these 
results are reliable. 

Since the duration for atmospheric stationarity is relatively short even when 
quantization techniques are used, the local information must be enhanced. Enhancement 
is obtained by using the statistical error to identify the significant transport information. 
Such results are shown in Figure 41, where the covariance is plotted as a function of the 
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Figure 41. Crossed-beam and anemometer comparisons. 
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convection velocity and compared to the anemometer velocities taken at the point of 
beam intersection. To establish the reliability of these results, the probability of retrieving 
the transport information [equation (39)] is calculated. In this case the probability was 

57 percent. In general, the probability of transport retrieval ranged from 50 to 58 percent 
in the clear-air turbulence experiments, where the transit distance ranged from 0 to 
1 00 meters. The empirical results show that the transport probability decreases as the 
transit distance increases. These clear-air turbulence measurements were all broadband 
processes with bandwidth ratios ranging from 4.1 to 6.3 and with an average frequency 
ranging from 0.01 to 0.02 Hz. The eddy size was determined to range from 73 to 
102 meters. The significance of these parameters for the clear-air measurements will 
become apparent when they are compared with the smoke transport in the atmosphere. 

The smokestack emission experiments are similar to the clear-air turbulence 
experiments in that both are atmospheric experiments and, in addition, the clear-air 
turbulence experiments are the limiting case for the smoke plume. However, the smoke 
plume is a dominant source of optical modulations, at least in the near field, and the 
smoke plume can, like the cold jet, be treated as having a point source origin. These 
differences permit the baseline of the radiometers to be aligned parallel with the axis of 
the smoke plume (Fig. 42) rather than perpendicular to the flow, as is done in the 
clear-air turbulence experiments. As a result of this test configuration, the transit height 
of the local information being retrieved is a variable which is uniquely determined from 
the transit time of the convective transport for a given set of geometric parameters. The 
reason for the different test arrangement is that for clear-air turbulence, which occurs at 
all altitudes, we wanted to isolate the altitude at which the transport phenomena are 
being measured; whereas, in smoke-plume measurements, which occur at a unique height, 
we want the velocities of transport and the height of transport [12]. Outside of these 
transit height considerations, the statistical retrieval of the convection speed is the same 
for the smoke transport as for the clear-air turbulence transport. 

Figure 43 shows a spatial mapping of the accumulative covariance curves resulting 
from smoke transport. The mean convective velocity was found to be about 3.2 m/sec at 
an altitude of 299 meters. Since we obtained a spatial mapping of these covariance curves 
along the streamline, we can obtain a dispersion velocity. In this case, the average 
frequency of the first covariance curve according to the wave packet algorithm was 
0.0491 Hz and decreased to 0.029 Hz for the last covariance curve. The bandwidth ratio 
was approximately 1.9 to 3.6, which is clearly smaller than that obtained in clear-air 
turbulence measurements. The mean dispersion velocity was determined to be 0.6 m/sec 
which means a lateral mixing rate of 0.2102. The eddy size increases from 34 to 

58 meters over a transit distance of 148 meters. This point of observation was about 
360 meters from the point of emission. This information for the convective and dispersive 
transport was incorporated into the smoke plume picture shown in Figure 44. The 
dispersion cone shown in this figure is in agreement with the observed average rate of 
dispersion. 
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Figure 42. Crossed-beam test arrangement for fan system. 

Comparing the smoke plume dispersive transport parameters for the average 
frequency and the eddy size with the parameters obtained for clear-air turbulence, we 
infer that the smoke plume is approaching the clear-air parameters. In fact, a linear 
estimate says that, at about 700 meters downstream of the source, the smoke plume 
should obtain comparable scales with the atmospheric turbulence. 

Our transport probability has an additional function in these measurements, since 
it serves as a check to verify that the data are from a common streamline. The transport 
probability ranged between 64 percent for the shortest transit distance to 51 percent for 
the maximum transit distance. These data show that the probability of retrieving 
transport information always decreased as the transit distance increased. It is interesting 
to note that the probability of retrieving transport information is higher for intersecting 
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Figure 44. Remote mapping of high stack emissions. 


beams in smoke plumes than in clear-air turbulence. This is logical since the smoke 
signatures are more intense. This higher probability means that the integration times 
required for correlating the smoke data can be shorter than those required for clear-air 
turbulence. 
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SECTION V. SUMMARY AND CONCLUSIONS 


The present investigation is the first study of dispersion and convection patterns 
of fluid transport phenomena using the crossed-beam correlation technique for the 
analysis of both stationary and nohstatiohary fluids. This new analysis scheme, as has 
been demonstrated, permits the separation of the dispersive parameters without extending 
significantly the data retrieval requirements. This separation of fluid transport was 
achieved by associating changes in the time domain with convection and changes in the 
frequency domain with dispersion. 

This new information retrieval scheme can be summarized as follows: 

1. The data source is the crossed-beam data collection system, which collects 
the temporal history of the optical fluctuations occurring in the field of view from the 
beams where they are in sufficient number to ensure spatial mapping. While the retrieval 
of the convection velocity does not require a spatial mapping, this mapping is highly 
desirable because it improves the statistical reliability of the results. On the other hand, it 
is mandatory that a minimum two-point mapping exist for the retrieval of the dispersive 
parameters. This spatial mapping is the only additional data collection requirement 
imposed by this analysis. 

2. This temporal history of the optical fluctuation is compressed into a 
correlation curve, from which the convection velocity and scales are directly obtained. To 
establish the stability levels and to employ statistical discrimination, it is necessary to 
subdivide the data record into pieces and to use the piecewise accumulative covariance 
technique. If the data tend to be of a nonstationary nature — that is, if short data 
recorders are used and thus rapid convergence is required — then it is advantageous to 
obtain piecewise independence by using the quantized accumulative covariance technique. 
This technique prevents any segment of the short data record from dominating the 
accumulative result. The stability tests were extended by the new concept of a transport 
probability, which affords an integral estimate for the minimum acceptable probability 
level for the retrieval of significant transport information and a check to determine 
whether the information is from a common streamline. 

3. The second-order cross-covariance is then analyzed for gradients in the 
average frequency composition. The gradients are the dispersive signature. This analysis is 
performed by suppressing the temporal information associated with the convective 
transport, which is the equivalent of changing from a fixed to a moving reference frame. 
In this way, the second-order covariance curve is transformed into a fourth-order 
covariance curve. The fourth-order covariance curve is curve-fitted, using regression 
analysis with the wave packet algorithm to obtain the frequency content in terms of the 
average frequency, the bandwidth, and the number density. This wave description is then 
used to obtain the dispersive transport, the lateral mixing scale, and the mixing rate 
characterizing the transport process. 
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4. The ratio of the bandwidth to the average frequency affords a systematic 
procedure for the classification of the various transport processes in accordance with their 
bandwidth ratio. 

This new information retrieval scheme was successfully applied to the analysis of 
transport phenomena in the glow discharge, cold jet plumes, and air quality 
investigations, where it was demonstrated that this procedure could meaningfully retrieve 
the dispersive and convective characteristic of the transport process. The most significant 
contributions of this analysis procedure are: 

1. The dispersion and convection phenomena were precisely defined from the 
equations of motion in terms of the normal and parallel fluxes of the stream tubes in the 
average velocity field. 

2. In the glow discharge, the ability to associate the longitudinal length scale of 
a moving stria tion to the length of the positive column, and the longitudinal length scale 
of a standing stria tion with the distance between striations was achieved. This length scale 
was also an index to the equilibrium of the discharge. 

3. The dispersion rate and the lateral mixing scales of the plume in the cold jet 
were predicted analytically. It was shown that this was in accordance with the empirical 
observations. 

4. In clear-air turbulence measurements, the procedures of suppression of 
nonstationarity and of peak identification were demonstrated, which resulted in the 
enhancement of the covariance curve. 

5. In the smoke plume studies, the convection and dispersion transport 
information was again retrieved from the crossed-beam data and was in agreement with 
the experimental conditions. 

From these demonstrations, we can conclude that this information retrieval 
procedure has successfully made an original contribution to the recognition and 
interpretation of dispersive and convective transport phenomena in both stationary and 
nonstationary fluid transport processes. 
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APPENDIX 


THE ENERGY DENSITY WAVE NUMBER SPECTRUM 


The objective of this appendix is to establish the requirements for a conservative 
wave number spectrum, which will be referred to as the energy density wave number 
spectrum. 

It can be shown [33] that the orthogonal form for the Fourier amplitude 
coefficients about the mean is 


m 


.a = 2 2 i 

^ i=0 


/ . 27rki . 27Tla 

i m 6. , , cos + 0 . sin 

y j, odd n j, even n 


m = n - 1 


(A.l) 


j 1, 2, 3, . . . , m 


‘I* i % odd 


whefe i is the time delay index, j is the amplitude index, k is the wave number, 
and 5j ^xxx Kronecker delta used to switch between the sinusoidal functions. 

The Fourier transpose from the frequency domain to the time domain is 


m 

J=1 


= / 27 rki 

a 6. cos 

m y j, odd n 


, . 27rki 

+ 6. sin 

J, even n 


(A.2) 


where the brackets indicate the frequency transpose of the covariance coefficient. The' 
last two. equations form a Fourier transform pair, which permits the interchange of 
information between the time and frequency domains. 

The odd presubscripted amplitude coefficients are commonly referred to as the 
cospectra (real part). 


, C — . a 

k m j odd m 


(A.3) 
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while the even presubscripted amplitude coefficients are referred to as the quadspectra 
(imaginary part), 


,Q = 

k m 


a 

even m 


(A.4) 


These amplitudes are normally converted to a polar coordinate system whose radius 
vector, the gain, is 


1.G = 

k m 


= 2 =2 

C j Q 

km km 


(A.5) 


and the phase angle is 


k^m 


arctan 


k^m 

k^m 


(A.6) 


whose function is similar to that of the time delay in the time domain. 

The restrictions imposed on the frequency description of a finite, discrete 
correlation curve by the orthogonal wave number representation are a bandwidth and a 
number, density. Specifically, the lower frequency limit is set by the total time interval 
subject to the transformation, while the upper frequency limit is set by the interval of 
time, At, between the covariance coefficients in accordance with the Nyquist 
uncertainty principle. The number of covariance coefficients, n , that describe the 
correlation curve sets the number of parameters in the frequency description. Since one 
parameter is required for the mean value of the correlation curve, there must be less than 
n-1 equally spaced amplitude coefficients. 


Next, we shall establish suitable requirements for the spectral normalization. To 
achieve this, consider the significance of the variance of the covariance 


curve. 



, in terms of the amplitude coefficients 



of the Fourier series. 


Let us obs erve that the variance under consideration is not the data 
variance, (Ax)^ [equation (19)] , commonly used in spectral normalization [8, 34] , but is 
the variance of the correlation curve that is actually being transformed into the frequency 
domain : 
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In accordance with equation (A. 2), the quantity in parentheses can be expressed as 
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which is the Fourier transpose. The above result is now substituted into equation (A.7), 
and the resulting matrix is reduced. When the arguments of the sinusoidal products are 
different, the uncorrelated terms are zero; thus, we obtain a diagonal matrix of the form 
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where the summations have been interchanged. 

Since a squared sinusoidal function over an integral number of waves is one-half, 
it follows that equation (A.9) reduces to the form 



(A. 10) 
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Obviously, the same result is obtained when the time lag range is extended to include 
negative time delays. Equation (A. 10) enables us to define an orthonormal wave number 
spectrum : 


1 




(A.ll) 


We can now define the terms describing the spectrum. When the amplitude 
coefficients are squared, they will be called intensity coefficients (same terminology as in 
geometric optics). Since the intensity coefficients describe the relative energy content in 
the covariance curve as a function of wave number, this spectrum is referred to formally 
as the energy density wave number spectrum. This is a conservative spectrum. 

The maximum number of linearly independent wave numbers in the 
representation for the frequency domain is (n-l)/2, where n is the number of 
covariance coefficients being transformed. This can be demonstrated by the following 
argument. The amplitude of the gain is 
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n + 1 
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j=0 •' 




27rkj ) 
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k = 1, 2, 3, . . . , n - 1 , 


(A. 12) 


where i stands for the square root of minus one. In complex notation, the amplitude of 
the gain can be expressed in terms of the cospectra and quad-spectra as 


= feC * 1 . (A.13) 

The amplitude of the gain [equation (A. 12)] can then be expressed in trigonometric form 
as 
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n + 1 
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27rkj 
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(A. 14) 


84 



For manipulative ease, we now write this equation [equation (A. 14)] in terms of the real 
and imaginary parts: 


k 
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1 

n + 1 


Yi .r cos 
j=0 ^ 



(A. 15) 
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(A. 16) 


Observe that when the wave number is zero (k = 0) , then the amplitude gain is 


G = 

k=0 n + 1 


s 

j=0 


(A. 17) 


which is the mean. This is why the wave number spectrum starts at one and has 
only n - 1 terms; i.e., one parameter of the spectrum is used for the mean value of the 
covariance curve. 

Next, consider the cospectra [equation (A. 15)] . The cospectra for the nth wave 
number is 
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(A. 18) 


Using the trigonometric angle-sum relation, the above equation can be written 


C = 
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(A. 19) 


where 


cos 27Tj = 1 


(A.20) 
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and 


sin 27rj = 0 


Substituting the last result into equation (A. 19), we obtain 




r cos 


27T 1 i 

n + 1 


(A.21) 


which is simply 


T C = , .C . 
k=n k=l 


(A.22) 


It logically follows then that, in general. 


k^ n-k+1^ 


1 ^ k ^ - 
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(A.23) 


The quadspectra [equation (A. 1 6)] has the same interrelationship as the 
cospectra. This is readily seen upon examination of the nth term that results from the 
application of the trigonometric angle-sum relation: 
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Again using the relations given in equation (A. 20), we obtain the fact that 


Q = Q 
k=n^ k=l^ 


(A.25) 
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or in general that 



n-k+1 


Q 


1 = k ^ ^ 
2 


(A.26) 


This derivation just proves the hmits placed on the wave number representation. 

Incorporating the normalization and limits just derived, we can write the Fourier 
transform for the wave number spectrum for the ampUtude density and intensity density 
representations as: 

1. Amplitude density representation, 
a. Cospectral amplitude density: 



(A.27) 


b. Quadspectral amplitude density: 



2 

(2n+ 



(A.28) 


where k = 1, 2, 3, ... , n/2. 



The lower case spectral coefficients indicate the normalization of the coefficients; 
whereas, the capital spectral coefficients are the unnormalized coefficients. The 
rectangular coordinate system discussed above has the following counterpart in the polar 
coordinate system; 

c. Gain amplitude density: 


k^m 





(A.29) 


d. The phase: 
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2. Energy density representation. 

a. Cospectral intensity density: 

V — 

kV km 

b. Quadspectral intensity density: 
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c. Gain intensity density (in polar form); 
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The energy density gain has the property that 


k=l 


It should be recognized that the averaging notation and the post subscript are carried 
over from the covariance coefficient to indicate the way this coefficient was obtained. 

These formulations for the spectra were tested by transforming the correlation 
curve in Figure 12 into the frequency domain. The results are shown in Tables A- 1 
through A-4 (the details of the calculations will be discussed later). From an inspection 
of these results, we can see that the majority of the calculations are insignificant. The 
next objective is to reduce the calculations of these insignificant coefficients. 

By using the properties inherent to the energy -density representation, a criterion 
for the minimum level of acceptable spectral intensity can now be developed. The 
spectral energy of the correlation curve in Figure 12 was concentrated in 9 of the 400 
wave numbers in this description, while the remaining 391 wave numbers contained less 
than 1 percent of the energy associated with the covariance. From a comparison of the 
known input parameters with these results, it is obvious that these 391 wave numbers 
represent the intensity level of calculation noise affording no insight into the behavior of 
the phenomena; therefore, a set of criteria to discriminate the significant energy bearing 
wave numbers from the noise bearing wave number is essential. 

The term “noise” is the key to the establishment of this cutoff level. A white 
noise process is a purely random one characterized by a time series, consisting of random 
impulses, that results in a cross-covariance curve of zero magnitude at all time lags and an 
autocovariance curve with a delta function at zero time delay. Since a preferred 
frequency would not result in this kind of correlation curve, the intensity of the gain 
must therefore be of the same magnitude for all wave numbers. Because the energy 
density representation has the property that the sum of the gain intensities is one 
[equation (A.35)] , it then, intuitively, follows that the intensity of the white noise gain 
is 


y 

white noise 


1 

n 


(A.36) 


This result represents the white noise level which is composed of an equal 
distribution of energy between all wave numbers; therefore, an energy less than the white 
noise level can be viewed as a lack of significant energy. Thus, the intensities that are less 
than the white noise level can be legitimately rejected as being below the intelhgence 
cutoff. 
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3.3-COCll 

0.0-30 ••-2.00c *• 

99.892> 

118 


14.75:‘03 •• 

:-.3j-3. 11 


:.C3. •• 

-3.-.3C-3;-: ♦ 

:.:o3 *• 

O-Oi’-Oll 

0.Q3G **-C.20: *• 

99.8923 

119 


14. 87530 •• 

-3.0(33^11 


}.co: •* 

0.--0C0:: * 

- -.233 •• 

j.OICOIl 

o.coo •*-r».coo ** 

99.8923 

12Ci 


la.OOOOO •• 

o.oooeu 


'.y-3. •• 

-C.'OOO'.- ♦ 

-.333 •• 

J.3C0-311 

0.3C3 **-0.03C *• 

99.8923 

121 


la. 125-3: •• 

-o.-:oo. ID 


:.-o- •• 

C* . lOOO • •; • 

c.':c3 •«- 

3.01.21 ' 

D.-303 **-:.2-2: ** 

99.8923 

122 


la. 26033 •• 

0.030M0 


O.OG- •• 

— o.-ooor-o • 

0.0C3 •• 

3.00-3--13 

3.*)3> •*-2,0y-2 •• 

. 99.8923 

12i 


15.376:C •* 

-o.cooao 


• ‘-.O-- •• 

o.cooo:*. • 

'J.-3-3 •• 

0.32-301 : 

0.3-3 **-c.oo: •• 

99,6923 

124 


15.6JJ3.3 •• 

-.000 -i: 


' • wOj •• 

-3.30001C « 

•’.JC3 •• 

3.3:001. 

0.030 ••-0.02: ** 

94.6923 

125 


15.6253. •• 

-C.OO-L'IO 


C.003 •• 

-.cococo 9 

3.303 •• 

0.3-3301.) 

O.'JOO **-;.33‘ •* 

•V9.8923 

12b 


15.76;:: •• 

0.300309 


— '3 L • • 

-3.0000 ’a • 

3.303 •• 

,.333009 

0.0-33 ••-!:. 33C •• 

99.8923 

127 


16.87533 •• 

-j, 033. 39 


l.co: •• 

2.COC3-1) • 

3.003 •• 

3.33';339 

•3.300 **-;.3-:: ** 

99.6923 

12e 


i6.:'0'jC3 •« 

C.0'3OC09 


‘.O-O: •• 

-O.f-0003- • 

0.-303 •• 

D.0C30-39 

0.000 ••-2,020 •• 

99.8923 

129 


16.12a:.’ V 

-o.ooo: J9 


:.,c-oc •« 

oioocoo; • 

C-.'5-33 •* 

0.0:-3339 

0,2)2 ••-:.301 •• 

79.6923 

153 


lb . 2a *• 

•3.000039 


-5.00- •• 

--.<•00330 • 

:.r-3 •• 

3.3--3039 

0.-2CO •*-:-. oo:- *• 

•79.8923 

131 


16.375.0 •• 

-C.003':09 


-.03. •• 

o.coooc: • 

0.0*33 •* 

0.00-3C09 

0.300 **-0.000 *• 

99.892) 

132 


lb. S3 j3v 4* 

o.3Co:L-a 


5.-3. •• 

-o.ioco:- • 

:.303 •• 

0.0-0C38 

0.00-) ••-0.J07 •• 

99.8923 

133 


l6.b253T •• 

-3.-30C: 38 


:.u-3C •• 

o.oooo:: • 

0.003 •• 

3.0-3208 

3. .30: **-1'.23:- *• 

99.8923 

134 


16.75::3 •• 

0.000.08 


.00. •• 

-o.cooo:: • 

-.333 *• 

3.0:0008 

0.0v>'3 **-0.020 •• 

97.6923 

135 


16.875:-.' •• 

-c.ooc:o8 


V.03-. •• 

-3. coco:;- • 

-;.303 •• 

3.010008 

0.23': 300 •• 

79.6923 

136 



C-.-3i33:08 


3.00- •• 

-0.---3C0O • 

C.*33D •• 

3.030038 

2.-300 ••-;). 3:: *• 

99.6923 

137 


17.126':-./ •• 

-:.300*'jb 


:.co- •• 

0.0000-?- * 

C.0C3 •• 

3.3f;-3D08 

0.000 ••-0.0-3C- *• 

•>9.6923 

138 


I7.26o03 •• 

0.330-;58 


:.oo.- 

-C.CV0053 • 

. 1* r, *• 4 9 

3.0. 3006 

3.032 •*-:,J2. *• 

99.6923 

139 


17.37a33 •• 

-3.000007 


•-.coo •• 

3.-.'303-3 • 

3.333 •• 

0.333CJ7 

o.c-3 •*-':,3&': •• 

99.6923 

143 


17.5o3u: •* 

0.000007 


<1,00- •• 

•O.C--00'3 • 

0.3C3 •• 

C.3C3337 

•-,3;o •*-•;. ouo •• 

97.6923 

141 


17.626:: «• 

-3.-300:07 


C.OOL- •• 

3.:'300- ; • 

-. 323 •• 

3.0;:037 

2.30-2 ••-3.03: •* 

y-).6923 

142 


i7.75::'w *• 

0. 000337 


C.OOi. •• 

-0.0000 - • 

C.333 •• 

0,3:-3037 

D.-320 ••-.-.30: •• 

99.8923 

143 


l7.o750y •• 

-3.0C0037 


*.00- •• 

3.2300-- > • 

0'.:D3 •• 

0.003007 

0.300 **-•:. 00: •• 

99.89^3 

144 


Id.'.Oj.'O •• 

0.000 :;'7 


->.o:- •• 

•0.-330': • 

c,:-3 •• 

0.0C3C37 

0.330 **••-. 302 •• 

79.8923 

145 


18.1 2 aoO •• 

-0. J003C7 


C.G30 •• 

0.0000”3 • 

-3.003 •• 

3.0CC007 

3.330 ••-0.:-3C’ •• 

99.8923 ■ 

140 


i8.25’;v: •• 

0.330-37 


O.C'OO •• 

-O.jOOO-'O • 

3.-3C3 •• 

.0,330007 

0.030 **--.23: •• 

79.8923 

147 


18.3750: 

-0.000:36 


•• 

C.)0C0*'-5 • 

3. •-OS •• 

3.0-30036 

3. -300 **-2.32C •• 

>•7.8923 

14B 


18.5C3C0 «« 

■3.000C06 


O.OOo •• 

-0.000030 • 

3.003 «• 

0.0-30006 

3.003 ••--.33-:- *• 

97.8923 

14V 


10.625C-3 «« 

•0.0000. 6 


C.C-- •• 

0.-3-3C-3*; • 

:.003 «• 

3.332006 

0.020 «*-3.30-: •• 

•79.8923 

153 


IB.76JOO *• 

:. 000336 


'/•o--. •• 

-o.-.cooo: • 

:.v33 •• 

3.020206 

2.030 •*-:.oo: •• 

99.8923 

151 


I8.o7a0i •• 

-C. 000006 


C.03- «• 

o.:ooo:o • 

0.303 •• 

0.000006 

o.coo ••-o.od:- *• 

9'y.8923 

152 


1-7., Jo--: •• 

3.300-*.J6 


-l.CO, •• 

-3..-303C3 • 

0.303 •• 

0.033C06 

0.000 ••-o.OUC *■ 

99.8923 

153 


19.12aL-; «* 

-3.3000-36 


o.c:: •• 

o.:oco-: • 

0.'-03 •• 

3.3:-J036 

0.003 **-:,30C •• 

79.8923 

154 


i9.26r-:o ** 

0.-3003:6 


".OOw •• 

-o.oocoro • 

3.303 •• 

0.303006 

0.000 **-3.000 *• 

97.8923 

155 


I7.375'.v •• 

-3.y003-:6 


-J.OO.- •• 

o.:ooo': • 

3.3-33 *• 

O.OvOCOb 

3.020 ••-3.32' •• 

99.8923 

156 


19.60::3 •• 

0.303006 


:-.co'- •• 

-0.330333 • 

3.303 •• 

3.3>3O06 

0.00-2 ••-0.002 

97.D923 

157 


19.b26.o «• 

-3.n00C05 


O.OOv •• 

c..-:300'Vj • 

3.003 •• 

0.0D3005 

O.'OO **-0.330 •• 

97.8923 

158 


19.75}:0 •» 

0.3000-5 


:.lo: *• 

-3.:-JO0l: • 

0,303 •• 

3.0:2035 

0.000 ••-2.')3'- *• 

99.8923 

159 


19.876/0 •• 

-0.-00035 


3.0-3- •• 

o.-:3oo-, / • 

-.003 •• 

o.o:oco5 

3.--30 ••-C.2-3 •• 

97.8921 

160 


2J.0Cj:;- *• 

o.03o::-5 


f.oo: *• 

-O.-iOCOv; • 

3.303 •• 

0.00D005 

O.C-3-3 **-3.302 •• 

99.6923 

161 


2j.l2a‘;: •• 

-•-.ooocos 


■ ‘.go: •• 

c./jOO:: • 

•3.333 •• 

0.000035 

C.C33 ••-U.333 •• 

99.8923 

162 


23.25JJ: •• 

:.3C0D;5 


i.yy: •• 

-3.;?co*7 • 

0.003 •• 

0.000035 

O.C3;3 **-3.2:3 ** 

99.8923 

163 


2J.37aCD •• 

-0.0000-5 


c.ooo •• 

3. coco-: • 

3. 303 •• 

0,330095 

3.000 **-^.001 •* 

99.8923 

164 


2-. 50-.: •• 

v.'j003-:/5 


1 .:-:; •• 

-0.:303'3 • 

•:*.oo3 •• 

3.320035 

0.333 ••-C-.2U-3 *• 

97.8923 

165 


20.62603 •• 

-C-. J0C.J36 


3.03- •• 

3.:oco. : • 

0.333 •• 

3.033005 

3.200 ••-C-.303 *• 

99.8923 

166 


23.7a3L0 •* 

3.030C -5 


«).o:: «« 

-c-.oc-oo-'i; ♦ 

3.003 •• 

3.030005 

3.000 ••-U.20C •• 

99.6923 

167 


23. 876-0 •• 

-0. JCO';35 


• ',.03. •• 

:'.::oo;- • 

0.003 •• 

3.3.-0035 

0.202 **-:.33'J *• 

97.8923 

16o 


2l.-0J':0 

0.000005 


2.033 •• 

-3,1000'.- • 

3.303 •• 

0.033035 

0.030 **—'3.031 •• 

79.8923 

169 


2i.U5j, •* 

-3. j330;-4 


*.OJw •• 

o.co-jo.*: • 

•3.003 •• 

0.023024 

0.-3O0 03-1 •• 

79.8923 

173 


21.26j-3 •• 

0.3LO"/;'4 


c.oo- •• 

-3..’ocO'-: • 

0.003 •• 

3.333034 

0.030 **-v‘.02-: •• 

99.8923 

171 


2l.i750‘3 •• 

-0.000034 


0.30- •• 

0. -03313 • 

•3.003 •• 

•3.3:0C34 

0.003 ••-0.300 •• 

99.8923 

172 


21.50J0J •• 

U.CJU'.rJ4 


O.OOJ •• 

-O..-300-’O • 

0.003 •• 

0,32DO-3* 

0.000 ••-2.03. •• 

97.6923 

173 


21.b27:3 •• 

-C.300:-34 


c.ooo •• 

C . ’Oi/O -0 • 

0.303 •• 

3.030004 

3.003 ••-:.o-: •• 

9-7.8923 

174 


21. 75C/LO *« 

3.00CC-34 


0.000 •• 

-3.3vOO'*>3 • 

0.303 •• 

3.0.0004 

3.000 **-:.30C- ** 

97.8923 

176 


21.87603 •• 

-0.0030.4 


C.30-. •• 

O.OOOO.j • 

0.303 •• 

3.033004 

0.2-30 ••—•). 00c •* 

77 . 8 92 3 

17o 


22..o:j-j •• 

3. 000304 


c.o-:: •• 

-3.’0C3:C • 

C-.3C3 •• 

o.o:ooo4 

0.-300 •*-j,20: •• 

-9.S923 

177 


22.1250C •• 

-3.0003-34 


O.OOw •• 

0.000030 • 

0.003 •• 

0.020034 

0.33'J **-:.0C3 *• 

99.8923 

178 


22.25*103 •• 

C.3C3--4 


V.030 •• 

-0. >>00' 3 • 

O.003 •• 

•3.000024 

0.033 **-O.JOC- *• 

97.8923 

179 


22.375-3 •• 

-3.000C--4 


c.oo: •• 

•3.00C0-':v • 

3.303 •• 

0.:33034 

3.-303 ••-2,03- •• 

97.8923 

180 


22.50-36 •• 

0.003334 


■;.ooo •• 

-C,1003'''3 • 

3.233 •• 

0.02D004 

0,03- *•-'-.-022 •• 

97.8923 

181 


22.62673 •• 

-0.-330-04 


c-.oo: •• 

3,:ooo:- • 

3.303 •• 

0.3C30C4 

0.030 **-0.03: •• 

99.8923 

182 


22.76JC'C •• 

j.C00^--;4 


7.03J •• 

-3..«J'-3--- • 

3.-33 •• 

3.02:0-4 

0.-3C0 •*-..J0-» •• 

99.8923 

183 

• . 

22.87503 •• 

-3.000 .■04 


c-.oc-: •• 

0. COCO' 3 • 

0.CO3 •• 

0. 000004 

0.003 **-u.2J: •• 

77.8923 

184 


23.303:- •• 

C.000*...4 


0.00- •• 

-o.-:ooo:: • 

0.003 •• 

3.300034 

0.003 ••-2.C0- *• 

79.6923 

185 


23.1253J •• 

-o.ooo:o3 


j.OOL- •• 

3.'33C0'-3 * 

o.:o3 •• 

0.31-2033 

0.030 **-3.2-: •• 

79.8923 

186 


23.2>j:0 •* 

C. 300-13 


c-.oco •• 

-3.300330 • 

C.-303 •• 

0.0-0003 

3.000 ••-O.OvO •* 

97.6923 

187 


23.376.3 •« 

-O.O30C33 


<3. 30- ** 

3.:-oco'.- • 

0.0-C3 •• 

3.0C3003 

0.-330 ••-),02: •* 

97.8923 

188 


23.50300 •• 

0.00'j0j3 


0.03- •• 

-3.00003- • 

0.033 •• 

3.3--3003 

3.000 ••--,00: ** 

99.8923 

189 


23.62aC-3 •• 

-3.330C-C3 


O.OJj •• 

0. -C?00v • 

:-.oo3 •• 

0.0DD003 

3. COO **-0.502 *• 

99.6923 

190 


23.75COO •• 

W.003333 


3. 03- •• 

-3.30303-J • 

C.0C3 •* 

0.000003 

- 0.000 **-2.00 * •• 

77.8923 

191 


23.876:3 *• 

-3.-30J- -33 


•:.oo- •• 

3.300:- • 

3.3C3 •• 

2.320033 

3.03-3 **-3.33; •• 

99.6923 

192 


24.00303 •• 

O.OOC-33 


C.02'- •• 

-c-.!-033:: • 

0.003 •• 

3.020033 

3,000 **-J.00'2 •• 

99.8923 

193 


24.126C0 •• 

-0.003-::3 


c.ooo •• 

3..-230-': • 

0.003 •• 

0.0w3033 

0.000 **-2.033 •• 

97.6923 

19a 


24.25300 •• 

&.300503 


C-.030 •• 

-0.000033 • 

0.003 •• 

0.0CO093 

9.030 **-2.000 *• 

99.6923 

195 


24.37500 •• 

-0.000- i3 


0.030 •• 

J.C-OOO'.J • 

0.003 •• 

0.030003 

O.COO **-2. 003 •• 

99.8923 

196 


24.50000 •• 

C. 0001-3 


V .'JOO •• 

--.V3031- • 

U.-C-3 •• 

0.0L0C33 

•3.000 **-i .:o> •• 

99.8923 

197 


24.625'-0 •• 

-C.003'.o3 


3.000 •• 

3.-0C033 • 

0.003 •• 

3.030033 


99.6923 

198 


24.75000 •• 

0.003J-33 


C.OOt- •• 

-3..00030 • 

0.020 •• 

3.00CC03 

0.000 *«-2.C2C- «« 

99.6923 

199 


24.87630 •• 

-C.000JC3 


V.OO- •• 

3.C30O-- • 

3.033 •• 

0.000003 

9.000 **-0.20. •• 

99. 

233 


25.--03CU •• 

o.o-'0:-3 


O.OOu •• 

-o.-300o:o • 

C.3C0 •• 

0.0COO93 

0.000 ••-' 



lcSJLII.<L# UUThuT iPECfKLM 


-fRc3UESCY RANGE » M^ TO 1,500 HZ 

NUMBER OF LwUAL AMPLITUDE COSINE WAVES > 21 
FRbUJENCY I..TERVAL 8ETHEEN WAVES - J.05 HZ 


MINIMUM ACCEPTAbLE GAIt'l INTENSITY 


Za.lZi-i 

< B . b 2 :> K j 

23.75. -,0 
zo.bts-j:* 

Zb.vc;-.:. 
2o.l2‘>.j; 
2o.25Jv-> 
26.27-j>- 
2b. aJ.-.- 
2o.u2ai-' 
2b* 7aw j 
2b.(37a 
27. jijtO 
27. 12 j. . 
27.2a-- . 
27.27av-> 

27.5.. . . 
£ 7 . b 2 a ,'i 
27. Taj. j 
2 7 . <J 7 a j 
2d. jO.' . 
2G. Ua 

2 d . 2 5 a J 
23. >7a 
2 0 . a V . i 
2 o.o 2 a.-j 
2o. 7av-’, 
2o.d7> -J 

2 Y . . J . , 

29.l2a..-< 
29.25”..- 
2V.a7a j 

2 y. Pv :-i. 7 
2-7.b2jjJ 
29. 7aaC0 
2-J.e7aJv 

Ai.l25.-,/ 
3O.2S0 j 3 
3i.375-;j 
3 j.5j:.: 
3:;.62a-:j 
ij.7pj-:o 
3C-.0/P..7 
3i.0wy'jj 
31.125GJ 
31. 2aj.j 

31.37300 
31.av-.10 

31.025. *0 
3U TaJOj 

31. d73v) 
3^. : 

32. 125,0 
•32.25jL-'v 

32. 'j7p., 0 
32 . aU j 
32.023^0 

i 

32. C750;- 
3a. 

33. 12a, 0 
33.2aVjj 
33.375.;-. 
33.5i/,» 0 

3 j. 62a-.‘i> 
33.75/,j 

33. fa75.7j 
>4 ..jo.7c: 
3«». 12a: ,■ 
3*f . 2 pjv 3 
3Y. 175CO 
ab. aOjJj 
3*t . b2a j j 
3‘f. 7ajvC- 
3>.. 87530 

3a. 125 } j 
3a. 2530, 
35.375 ,,3 
3a. aJ jj-j 

35.025.. ' 
35.75jv 0 
33. 87aC: 
ib, 3 jvO J 
3b . i ii'Jv 3 

30.25. -J-: 
3o. J7ajj 
3o.aJor-: 
3o. 62a 3 j 

30.75. jO 
3b.B75j'j 
37. JO ‘.-J 
37. 1250J 
3 7.25 jC-b 

37.37300 
37.aOOjv 


COSPFCTRa • 

AMl'LiruuE «lUTt.lSlTY* 


0..00. j2 

-o.go-o;j2 
L.j?o:,2 
-j.'JOJ -2 
-,.OjO.J2 
-j. :oo, :Z 

-:.-o-'o-.j2 
:.j.- j'.2 
-0. J-.V0...2 
C.OOC. j2 
.-2 

V..-CC 2 

-0.000. j2 

-J.jjO •j 2 
-j.00-j.--2 
.i,OC-0-,2 
-O.C' J.:-2 
j.coo: 2 
-‘0.o:o;.-2 

C.C-'C -2 
-O.b.'O. J2 
O.-j-O -72 
-^>■.00 j2 
J00M2 
-.'.OCO .2 
(...OO, ;2 
-w.c:o :2 
O.bOO-. j2 

• :,>-oo..2 
-j!000..-2 

f-. /j-‘. .-I 
- j.OOO^ U 
.••.0-OJJ..I 

U.JOO'-Ol 
-C.OOO-.-l 
c.-ooD. :-i 

-C. jUCv >l 

:.ooo-'oi 

-j.-OOO .,-1 

o.ooo..;'i 
-:.jC-o. ; i 
: .c-ooc :i 

-O.JOOCOl 

o.ooo::i 

-O.OOO-.l 
0.000«*);l 
-C»b’jO J.-l 
O.vOOi .1 

-o.co.oc:-! 

O.OO.J.- jl 

- J.OJj ..I 
-J.'JOOCOl 

-j.OOO Ji 
O.OOO *;l 
-O.OOuC-Ol 
'j.OO-Ov.l-' 
. 700 rj I 

o.o-:oj:i 
-U.COOJjl 
j.'JOOO-l 
-C.OOOl Cl 
O.OOO.C-l 
-j. jOCC-l 
O.OOOCOl 
-O.OCOC-M 
-o.uoooji 
-c. jcc-:j1 
C.OOOIOI 
-C..JO-0-; ji 
•:.jo-o.-‘i 
-o.coo: :^i 

O.COOcjl 

- j,OOC'*-C-l 
O.OOO..JI 

-o.oooon 

O.OOOt-'.l 

-o.coo..'. 1 
0.000 <01 
-O.-JO-O^ 01 
0.000.' -01 
-C-.-OOOOJl 

o.ooocn 

-C.uOO.Ol 

j.c-Dor.«i 

-c.ooor.'i 

j.OOO.jjl 

-o.ooc-:n 

iO.OOOMl 

.-L 

c.-30o:-oo 

-c-.ooo:o:- 

O.OCOOjJ 


3. 0.30 ; 
-0..003':- 
0. ojj :• 
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The intelligence cutoff criteria, along with the conservative property associated 
with the energy density representation, are the ingredients for a spectral inteUigence 
bandwidth which prevents blind data ingurgitation. When the sum of the gain intensities 
is 


k 


E 

k=l 


cut-off 


kv 


1 


y 


white noise 


(A.37) 


then, by the conservative property [equation (A.35)] afforded by this representation, we 
know that the remaining gain intensities will not afford any additional significant 
intelUgence. This implies that we have found the upper limit of the intelligence 
bandwidth, and that the spectral calculations can be safely terminated without a possible 
loss of significant information. The lower limit of the intelligence bandwidth is set by 
where the spectral intensities first exceed the intelUgence cutoff. 

Applying the concepts developed in this appendix to the transformation of a 
correlation curve obtained from atmospheric turbulence data, we obtained the energy 
density spectrum shown in Figure A-1. Use of the traditional Fourier spectrum analysis 


INTENSITY, <y>mi PERCENT 



Figure A-1. Atmospheric energy density spectrum. 



techniques [8,34] would require an explanation of the entire spectrum. Use of the 
intelligence cutoff criteria isolates the spectrum into two regions: from 0.01 to 0.26 Hz 
and from 1.0 to 2.4 Hz. 

The significant point which has been established is that the intelligence cutoff 
criterion affords a systematic noise suppression which enhances the spectral intelhgence. 
The intelhgence termination point, on the other hand, is just a labor-saving scheme, 
permitting us to reduce the amount of wasted time and energy in searching regions in 
wave number space void of significant information. In this example, computational time 
could be reduced by one-third if we stop calculating the spectral coefficients after we 
reach the intelligence termination point. 


George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama, October 20, 1972 
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